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O | Abstract 

We first give an overview of the shell-correction method which was developed by V. M. 
Q | Strutinsky as a practicable and efficient approximation to the general selfconsistent theory of 

finite fermion systems suggested by A. B. Migdal and collaborators. Then we present in more 
detail a semiclassical theory of shell effects, also developed by Strutinsky following original 
ideas of M. Gutzwiller. We emphasize, in particular, the influence of orbit bifurcations on 
shell structure. We first give a short overview of semiclassical trace formulae, which connect 



the shell oscillations of a quantum system with a sum over periodic orbits of the corresponding 
classical system, in what is usually called the "periodic orbit theory". We then present a 
case study in which the gross features of a typical double-humped nuclear fission barrier, 



classical system, in what is usually called the "periodic orbit theory". We then present a 

m 
oo 

including the effects of mass asymmetry, can be obtained in terms of the shortest periodic 
orbits of a cavity model with realistic deformations relevant for nuclear fission. Next we 
investigate shell structures in a spheroidal cavity model which is integrable and allows for far- 
going analytical computation. We show, in particular, how period-doubling bifurcations are 

•i-H . 

^ ■ closely connected to the existence of the so-called "super deformed" energy minimum which 

h : 

corresponds to the fission isomer of actinide nuclei. Finally, we present a general class of radial 
power-law potentials which approximate well the shape of a Woods-Saxon potential in the 
bound region, give analytical trace formulae for it and discuss various limits (including the 
harmonic oscillator and the spherical box potentials). 
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1 INTRODUCTION 



This paper is devoted to the memory of A. B. Migdal. Our first aim here is to review the shell- 
correction method (SCM) which was first introduced by Strutinsky on a phenomenological basis 
[1] and then microscopically founded [2] on Migdal's theory for strongly interacting finite fermion 
systems [3]. Our second aim is the discussion of a semiclassical theory of shell effects, using the 
so-called periodic orbit theory (POT) (see [4] for an introductory text book). It provides us with a 
nice tool for answering, sometimes even analytically, some fundamental questions asked [5, 6, 7] by 
Strutinsky: Why are nuclei deformed? What are the physical origins of the double-humped fission 
barrier and, in particular, of the existence of the isomer minimum? His idea was to use the POT for 
a deeper understanding, based on classical pictures, of the origin of nuclear shell structure and its 
relation to a possible chaotic nature of the nucleons' dynamics. We shall present some applications 
of the POT to nuclear deformation energies and discuss in more detail the relation of bifurcations 
of periodic orbits with pronounced shell effects. 

According to the SCM, the oscillating part of the total energy of a finite fermion system, the so- 
called shell-correction energy SE, is associated with an inhomogeneity of the single-particle (s.p.) 
energy levels near the Fermi surface. Its existence in dense fermion systems is a basic point of 
Landau's quasi-particle theory of infinite Fermi liquids [8, 9], as extended to self-consistent finite 
fermion systems by Migdal and collaborators [3, 10]. This is schematically illustrated in Fig. 1, 
where the s.p. level spectrum of a bound nucleus is shown in two extremal situations. Depending 
on the level density at the Fermi energy - and with it the shell-correction energy SE - being a 
maximum or a minimum, the nucleus is particularly unstable or stable, respectively. This situation 
varies with particle numbers and deformations of the nucleus. In consequence, the shapes of stable 
nuclei depend strongly on particle numbers and deformations. This is illustrated in Fig. 2. Here 
the shell correction 5E of the neutrons is shown as function of the neutron number N and the 
deformation parameter rj of a Woods-Saxon potential [11] with spheroidal shape, i] being the ratio 
of the semi-axes. If we fix the neutron number N, e.g. N = 150, and increase the deformation rj, we 
meet the first minimum (ground state) at about rj ~ 1.25 and the next one (isomeric state) at much 
larger deformations rj ~ 1.9 — 2.1. The experimental data corresponding to these deformations are 
shown in Fig. 2 by the heavy dots. 

The SCM was successfully used to describe nuclear masses and deformation energies and, in 
particular, fission barriers of heavy nuclei. For an early review by Strutinsky's group, in which also 
the miscroscopic foundations of the SCM are discussed, see Ref. [12]. (We refer to Sect. 3.2 for a 
further discussion of fission barriers.) 



In Sect. 2, we will give a short review of the SCM and its foundation on the basis of a selfcon- 
sistent theory of finite interacting fermion systems. 

Sect. 3 is devoted to the semiclassical theory of shell effects. The POT is based on Gutzwiller's 
semiclassical trace formula for the level density for a Hamiltonian system with isolated orbits [13] 
and its extensions to systems with continuous symmetries [4, 5, 14, 15, 16, 17]. It allows one to 
relate both the oscillating part of the level density and the shell-correction energy of a quantum 
system to the shortest periodic orbits (POs) of the corresponding classical Hamiltonian system. 
Thus, one can often explain pronounced shell effects by the role of particular short POs. As an 
early example, taken from Ref. [5], the heavy bars in Fig. 2 are the predictions of the POT for the 
loci of the ground-state minima, using the shortest POs in a spheroidal cavity. Bifurcations of POs 
under the variation of a deformation parameter or the (Fermi) energy can have noticeable effects 
for the shell structure [5, 6, 7, 18, 19, 20, 21]. 

In Sect. 3.1, we will present the structure of semiclassical trace formula and discuss a general 
method of treating bifurcations in the POT, using the catastrophe theory of Fedoryuk and Maslov 
for caustic and turning-point problems [22, 23, 24, 25]. 

In Sect. 3.2, we review the semiclassical description [26] of a typical nuclear fission barrier in 
terms of the shortest periodic orbits, employing a cavity model with the realistic shape parame- 
terization developed in [12] . In particular, the effect of left-right asymmetric deformations on the 
height of the outer fission barrier will be discussed. Isochronous bifurcations of the shortest or- 
bits are treated here in a uniform approximation employing a suitable normal form for the action 
function. 

In Sect. 3.3, we use the spheroidal cavity [5, 19, 20, 27] as a simple integrable model that allows 
to study semiclassically the shell structure related to the 'super-deformed' energy minimum which 
in realistic actinide nuclei corresponds to the fission isomers. Strutinsky's prediction concerning 
the importance of the enhancement of the shell structure owing to period-doubling bifurcations of 
three-dimensional POs from simple equatorial (EQ) orbits in the spheroidal cavity model [5] will 
be discussed. 

In Sect. 3.4, we shall study a radial power-law potential [28], which is a good approximation to 
the familiar Woods-Saxon potential for nuclei in the spatial domain where the particles are bound. 
We shall establish generalised trace formulae for this potential and discuss various limits to other 
known potentials. 

The paper is summarized in Sect. 4, where we also present some conclusions and plans for future 
research. Some technical details of our POT calculations are given in the Appendix. 



2 THE SHELL-CORRECTION METHOD 



In 1966, Strutinsky achieved a far-reaching break-through [29], following basically Migdal's theory 
of finite fermion systems [2, 3]. Until then, many attempts had been made to incorporate quantum 
shell effects in the calculation of nuclear deformation energies. But they all failed in reproducing the 
fission barriers of actinide nuclei and details such as, e.g., the left-right asymmetry of the nascent 
fission fragments [30]. Summing the s.p. energies of a deformed shell model (like the Nilsson model 
[31]) up to the Fermi energy failed at larger deformations. There was a need to renormalize the 
wrong average part of the s.p. energy sum. Knowing that the smooth part of the nuclear binding 
energy could be well described by the phenomenological LDM [32, 33, 34] (or droplet model [35]), 
Strutinsky wrote the total nuclear energy as [1, 12] 

E to t — E LDM + SE , (1) 

where Eldm is the LDM energy and SE the so-called "shell- correction energy" which contains the 
fluctuating part of the s.p. energy sums for the neutrons: 

N/2 N/2 

SE = J2En-(J2 E n), (2) 

71=1 71=1 

and similarly for the protons. Both parts of the total energy (1) depend on the neutron and proton 
numbers N and Z, as well as on the nuclear deformation which has to be suitably parameterized both 
in the LDM and the shell model. While it had been a wide-spread belief that the shell-correction 
SE was important only for spherical nuclei and would vanish at larger deformations, Strutinsky was 
convinced that shell effects play an important role also at larger deformations and lead, in fact, to 
new magic numbers corresponding to deformed systems with increased local stability [1]. 

For the determination of the smooth parts (J2n=iEn) in (2), Strutinsky designed a very in- 
genious averaging method [1, 29] which has been termed the "Strutinsky averaging (or smoothing) 
method". It consists of a Gaussian convolution of the s.p. energy spectrum, modified by the so-called 
"curvature corrections" in such a way that the result does not depend on the energy averaging width 
7 (at least within a finite interval of 7 of the order of the main shell spacing) and at the same time 
reproduces the true average level density which is found, e.g., from the extended Thomas- Fermi 
(ETF) model, or from the Weyl expansion in the case of billiard or cavity models (see, e.g., [4], 
Chap. 4, for details). 

Strutinsky applied the SCM to the calculation of fission barriers [1, 29], employing the Nilsson 
model. For a typical actinide nucleus he obtained a second minimum in the deformation energy, 
lying above the ground state by about 3 MeV, at a deformation much larger than that of the 



ground state. This was, in fact, the physical explanation of the fission isomer which had been 
known experimentally [36] since 1962, but not understood theoretically. Strutinsky presented this 
result at the Symposium "Nuclides far off the stability line" in Lysekil (Sweden) [37] in 1966 and, 
being a member of Migdal's theory group, he immediately became famous. He was then invited 
to the NBI in Copenhagen, in order to extend his calculations of fission barriers on a larger scale, 
employing more realistic shell-model potentials, which led to the team work published in Ref. [12]. 

The SCM was soon taken up by many groups of the international scientific community and 
often also called the "microscopic-macroscopic" method [35, 38]. Various combinations of nuclear 
shell models and liquid drop(let) models were used. Still today, the shell-correction method is being 
used world wide for calculations of nuclear masses and deformation energies, and it persists to yield 
the most accurate nuclear mass tables, ground-state and isomeric deformations and, in particular, 
fission barriers. 

The Strutinsky averaging used for the second term in (2) is done in energy space and leads, 
within a sufficient numerical accuracy, to the same results as the ETF model ([39], see also [4], Ch. 
4.7). However, it is formally not completely consistent with the particle-number averaging that is 
implicitly done in the standard least-square fits to the LDM which defines the smooth part of the 
total nuclear energy (1). Already early, alternative particle-number averagings were investigated 
in [40]. Using semiclassical POT arguments, the difference between energy and particle-number 
averaging was understood [41, 42] as a symmetry correction which becomes especially significant for 
spherical nuclear shapes (particularly in the harmonic oscillator model). Smaller discrepancies for 
the deformed Fermi systems in the resulting shell-correction energies persist, however, as discussed 
in [41, 42, 43]. This point is thus still an object of current debate. 

The decomposition of the total nuclear energy into a smooth and an oscillating part in (1) 
may at first glance look like a rather phenomenological ansatz. In particular, since the oscillating 
part is taken from the sum of occupied s.p. energies of the nucleons, one may argue that it cannot 
be correct, since this sum is well known not to represent the total energy in a self-consistent 
microscopic theory (where it double-counts the potential energy if the interaction does not depend 
on the density). However, it was soon realized that (1) is nevertheless correct even within a self- 
consistent microscopic theory. This was pointed out, amongst others, by Bethe [44] who therefore 
termed (1) the "Strutinsky energy theorem". It can be rigorously proved that the oscillating part 
of the correct total energy is, indeed, contained in the s.p. sums; the proof is simply based on the 
variational principle that governs the self-consistent mean- field theories (see, e.g., [4], App. A. 3), 
and it applies also to density-dependent nuclear interactions [44]. This was soon demonstrated by 



the Strutinsky group [2] to hold also within the Migdal theory [3]. At this point it may be worth 
mentioning that the Migdal theory has been used explicitly to calculate energy shell corrections for 
doubly magic nuclei, investigating the so-called "lead anomaly" (i.e., in particular for nuclei around 
208 Pb) [45]. 

The Strutinsky energy theorem (1) was later tested numerically using the Hartree-Fock (HF) 
approach with effective Skyrme interactions. By extracting a self-consistently Strutinsky-averaged 
part of the total HF energy (which is ideally represented by the LDM energy E LDM ), the first-order 
oscillating term was, indeed, found to be correctly represented by the shell-correction energy (2) 
evaluated in terms of the energy spectrum of the averaged self-consistent HF potentials (which are 
ideally represented by the shell- model potentials). The remaining higher-order terms were found to 
be relatively small, less than ~ 1 MeV in all (not too small) nuclei [46]. 

The shell-correction method is thus a well-established practical approximation to a self-consistent 
microscopic theory and applicable to any bound many-fermion system. It has, e.g., been applied 
to metal clusters by many groups (see [47, 48, 49, 50] for a few representative references). For 
semiconductor quantum dots, the Strutinsky energy theorem and shell-corrections to the Coulomb 
interaction energy were discussed in [51]. 

3 SEMICLASSICAL THEORY OF SHELL STRUCTURE 

3.1 SEMICLASSICAL TRACE FORMULAE 

In the semiclassical trace formula of Gutzwiller [13] that connects the quantum-mechanical density 
of states with a sum over POs of the classical system, divergences arise at critical deformations where 
bifurcations of POs occur or where symmetry breaking (or restoring) transitions take place. At these 
points the standard stationary-phase method (standard SPM or SSPM), used in the semiclassical 
(asymptotical) evaluation of the trace integrals, breaks down. Various ways of avoiding these 
divergences have been suggested using uniform approximations (see [4], Sec. 6). Presently, we shall 
discuss the evaluation of trace integrals in the phase-space representation, following mainly the 
presentations in [18, 19, 20, 21]. 

The essence of a semiclassical trace formula is to relate the oscillating part of the level density 
of a quantum Hamiltonian system to the periodic orbits of the corresponding classical system. The 
level density g(E) can be obtained from the semiclassical Green's function by taking the imaginary 



part of its trace in phase space variables, see [21], also the references therein, 
g(E) = J2 S (E-E n )^ReJ2j^f 3 nE-H(r",p" 

x | r | ^ exp 
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- <P CT — l — fi CT 



(3) 



where the delta function imposes the energy conservation of the particle motion along the classical 
trajectory (CT) in the potential well of the Hamiltonian H{r, p). $ CT is the action phase, 

®ct = ScrP, r", E) - p" ■ (r" - r') ; S CT (r>, r", E) = £ dr ■ p(r) . (4) 

The second equation yields the action along the CT. The Maslov phase ficr is determined by 
the number of caustic and turning points within the catastrophe theory by Fedoryuk and Maslov 
[21, 22, 23, 24]. In (3), J7 cr (p / _ L , p") is the Jacobian of the transformation from the momentum 
perpendicular to the CT p' ± at the initial point r' to the perpendicular momentum p" at the final 
point r" of the CT. One may divide the level density into a smooth and an oscillating part: 

g(E)=g(E) + 5g(E). (5) 

The smooth part g (E) is, to leading order in 1/h, given [52] by the direct (zero-time or zero-length) 
trajectories which yield the Thomas-Fermi (TF) results. Except for one-dimensional systems, higher 
h corrections contribute to it which may be calculated within the ETF model (or the Weyl expansion 
for billiards) [4], so that we can make the identification 

g(E)=g ETF (E). (6) 

For the calculation of the oscillating component 5g(E), we apply the stationary phase method 
to the integration over the phase space variables. The stationary phase condition reads: 



P") = , (7) 



and is nothing but a condition for CT to be periodic. In the presence of continuous symmetries, 
the stationary points form a family of periodic orbits (POs) which cover a (/C+l)-dimensional 
submanifold T po of phase space (denoted T po in Ref. [4] and further references cited therein), whereby 
K, is the degeneracy of the PO family. The integration over T po must be performed exactly. In 
integrable systems, it is advantageous to transform phase space variables from Cartesian to action- 
angle variables (see, e.g., [16, 17]). Then the action $ is a function only of the action variables, 
and the integrations over the cyclic angle variables are exactly carried out. Integrating over the 
remaining action variables using the standard SPM (under the existence of additional symmetries 



like SU(3) or 0(4), one of which can also be performed exactly), one obtains the so called Berry- 
Tabor trace formula [16] . 

For solving bifurcation problems in integrable and non-integrable systems, more exact integra- 
tions are required. In the SPM, after performing the exact integrations over T po , one uses an 
expansion of the action phase <3>ct m phase space variables £ = {r', p"}j_ perpendicular to T po in 
the integrand of (3) over £ near the stationary point £*, 

*CT(0 = % + l^oiCm - C) 2 + iKoiOti - C) 3 + • • • , %o = $ CT = *CT(0 • (8) 

To demonstrate the key point of our derivations of the trace formula, we consider here only one 
(one-dimensional) variable, called £ again, from the phase space integration variables in (3), on 
which we meet a bifurcation (catastrophe) point in applying the SPM. (We shall give comments if 
this might lead to a misunderstanding.) In the standard SPM, the above expansion is truncated at 
the 2nd order term and the integration over the variable £ is extended to ±oo. The integration can 
be performed analytically and yields a Fresnel integral, see e.g. [4], Sect. 2.7, Eq. (2.169). 

However, one meets singularities using the standard SPM which are related to zeros or infinities 
of $p (£*) (or of eigenvalues of the corresponding matrix in the case of several integration variables 
£) while $po(£*) remains finite in the simplest case. These singularities occur when a PO (isolated 
or family) undergoes a bifurcation at the stationary point £* under the variation of some parameter 
(e.g., energy or deformation). The Fresnel integrals of the standard SPM sketched above will then 
diverge. In order to avoid such singularities, we observe that the bifurcation problem is similar to 
the caustic singularity considered by Fedoryuk within the catastrophe theory [22, 24], adopted to 
its specific position at the edge of the phase-space volume accessible for the classical motion (see 
also Appendix A in [21]). Therefore, we employ what we call the "improved SPM", in short: ISPM 
[18, 19, 20, 21]. Hereby the integration over £ in (3) is restricted to the finite limits defined by the 
classically allowed phase space region through the energy- conserving delta function in the integrand 
of (3). The expansion (8) of the action phases and similarly of the amplitudes in (3) is generally 
used up to the second- and zero-order terms, respectively, and if necessary, to higher order terms 
in £ - £*■ 

In the simplest version of ISPM, the expansion of the phase is truncated at 2nd order, keeping 
the finite integration limits £_ and £ + given by the accessible region of the classical motion in (3). 
It will lead to a factor like 

e^° /h £ + exp [JLSJ, (£ _ d £ oc -J= e^°' h erf [Z_, Z + ] , (9) 



where erf(zi,z 2 ) is the generalized error function with complex arguments 



erf z 2 ) = -L r e~*dz, Z ± = (£ ± - f) J ~^% ■ (10) 



Note that the above expression (9) has no divergence at the bifurcation point where $p (£*) = 0, 
since the error function (10) also goes to zero linearly in \JW p r [cf. the second equation in (10)], 
which keeps the result finite. [For the case of several variables £ for which we find zeros or infinities 
of eigenvalues of the matrix $p (£*), we diagonalize this matrix and reduce the Fresnel-like integrals 
to products of error functions similar to (9).] 

This procedure is proved to be valid in the semiclassical limit by the Maslov-Fedoryuk theorem 
[22, 23, 24]. In this way, we can derive contributions from each periodic orbit free of divergences 
at any bifurcation point, and the oscillating part of the level density can be approximated by the 
following semiclassical trace formula: 

5g(E) ~ 5g scl (E) = J2 5 9po( E )> with ( n ) 



po 



6g£(E) = RelA po (E)exp 



1 

jS po (E) - io- po -n/2 - i(p d 



The sum is over all periodic orbits (isolated or families) (po) of the classical system (besides of the 
POs which belong to the same family on which the summation was already performed). S po (E) = 
§ p dr are their action integrals, the amplitude A po (E) (which in general is complex) is of the order 
of the phase space volume occupied by CT, and the factor given in Eq. (9) which depends on the 
degeneracies and stabilities of the POs, respectively. a po is called the Maslov index. 4> d is an extra 
phase that depends on the dimensionality of the system and degeneracy /C of the PO manifold. (4>d 
is zero when all orbits are isolated (/C=0), as defined in [13]). The sum in (11) is an asymptotic 
one, correct to leading order in 1/h 1 ^ 2 , and in non-integrable systems it is hampered by convergence 
problems [13]. For systems in which all orbits are isolated in phase space, Gutzwiller [13] expressed 
the amplitudes A po (E) (which are real in this case) explicitly in terms of the periods and stability 
matrices of the POs, see some examples below. His trace formula has become famous, in particular 
in connection with "quantum chaos" [13]. Notice that according to (9), any more exact integration 
in (3) over a bifurcation/catastrophe variable £ of the improved SPM leads to an enhancement of 
the amplitude A po in the form of a maximum in the transition interval from the bifurcation point 
to the region of the asymptotic (SSPM) behaviour of A po . The hight of this maximum is of order 
l/h 1 / 2 as compared to the result of the standard SPM integration (integrable or non-integrable 
system; see more specific examples in Sects. 3.3 and 3.4). 



The trace formula (11) thus relates the quantum oscillations in the level density to quantities 
that are determined purely by the classical system. Strutinsky, in his search for simple physical 
explanations of shell effects, realized that this kind of approach could help to understand the shell 
effects in terms of classical pictures. However, in the application to nuclear physics, Gutzwiller's 
expression for the amplitudes A po (E) could not be used, because they diverge when the POs are 
not isolated in phase space. This happens whenever a system has continuous (e.g., rotational) 
symmetries, and hence for most typical shell-model potentials (except in non-axially deformed 
situations). Gutzwiller's theory was extended to systems with continuous symmetries in [15] (see 
also Chap. 6 of [4] for details). 

Trace formulae for systems with all kinds of mixed symmetries, including the integrable cases, 
were also developed later by various other authors. The treatment of bifurcations is still an on-going 
subject of current research. Uniform approximations were constructed for orbit bifurcations and 
symmetry breaking under the variation of the energy or a potential parameter; references to most 
of these developments are given in [4], Chap. 6. 

Strutinsky has, however, not only the merit of extending Gutzwiller's approach to realistic shell- 
model potentials, but he and his collaborators also extended the semiclassical approach to the 
description of bound many-fermion systems in the mean-field approach. In [15] it was shown that 
for such systems the shell-correction energy 8E (2) (for one kind of particles) is given semiclassically 
by a similar-looking trace formula: 



The difference to the trace formula (11) for the level density is that here the amplitudes and actions 
are to be evaluated at the Fermi energy E F of the considered particles, and the appearance of an 
extra factor (h/t po ) 2 , where t po are the times periods of the particle motion along the PO with M 
repetitions, t po = MT po , T po being the primitive (M = 1) period of PO. This extra factor brings a 
natural convergence to the sum, different from that in (11): orbits with longer periods contribute 
less to the shell-correction energy (besides of amplitude enhancement of shell structure discussed 
above and further below). To obtain the energy shell correction 5E (12) as function of the particle 
number N, one can use the standard relation of the Fermi energy Ep to the given particle number 



The level density g(E) can be approximated here with Eqs. (5) and (6) using the trace formula Eq. 



5E ~ 5E scl = Re ^{h/t po ) 2 A po {E F ) exp -S po {E F ) 



ia po 7i/2 - i(j) d 



(12) 



po 



N, 




(13) 



(11) for its oscillating part. 1 

It is the beauty of this approach that gross-shell effects in 5E often can be explained semiclas- 
sically in terms of a few of the shortest POs (with smallest periods) in the system. The earliest 
application of this idea by Strutinsky's group was given in [5] , where the correct slopes of the stabil- 
ity valleys rj(N) of nuclei in a plot of 5E versus particle number N and deformation 77 was correctly 
reproduced using the shortest orbits in a spheroidally deformed cavity at small 77, as shown in Fig. 
2. A further example will be given in Sect. 3.2. However, if one wants to study finer shell struc- 
tures or some specific situations at large deformations, longer orbits have to be included. Hereby, 
bifurcations of POs can play a crucial role, as it will be exemplified in Sects. 3.3 and 3.4. 

For billiard systems, it is advantageous to express the level density as a function not of energy 
E but of the wave number k = y / 2mE/h 2 , where m is the mass of the particle. By virtue of the 
simple k dependence of the phase function S po /H in (11) for such systems through the product of 
wave number k and orbit length L po , the Fourier transform of the level density g{k) with respect to 
k yields directly its length spectrum and thus provides a nice tool to examine the quantum-classical 
correspondence. The Fourier transform of the semiclassical level density gives, indeed, 

F(L) = J g(k)e tkL dk ~ F (L) + ^ A po 5{L - L po ) (14) 

po 

which has peaks exactly at the lengths of the periodic orbits L = L po with heights proportional to 
their amplitudes A po . F (L) represents a Fourier transform of the smooth TF level density and has 
a peak at L — 0. Taking the Fourier transform of the exact quantum spectrum of the Hamiltonian 
thus reveals the classical periodic orbits, and the heights of the Fourier peaks give information about 
their semiclassical amplitudes. 

It is often instructive to study a coarse-grained level density which is defined by convoluting the 

exact sum of delta functions in (3) with a normalized Gaussian function of finite width T: 

1 r°° 1 
g T {E) = — / dE'g(E') exp[-(£ - E'f/^\ = —= £ exp[-(£ - E n f /r 2 ] . (15) 

Convoluting the semiclassical trace formula (11), one may expand the action S po (E') in the phases 
up to linear (and the amplitudes A po (E') up to zero) order terms in E' — E, as long as T is small 
with respect to the distance between the gross shells near the Fermi surface, see [15]. In this way, 
one obtains the Gaussian-averaged trace formula 



^r,sci(£) = S>S( E ) ex P |- (^) ] • ( 16 ) 



1 For practical purposes, we use the coarse-grained expression denned in Eq. (17) below with a small value of 7 to 
which Ep is quite insensitive. 



For billiards, a corresponding Gaussian averaging with a width 7 over the wave number k leads 
(with better accuracy owing to the linear dependence of the action on k) to the averaged trace 
formula 



where L po is the length of the PO, and the relation between the two smoothing widths is T ~ 2^E/k. 
We see that, depending on the smoothing width T (or 7), longer orbits are automatically suppressed 
in the above expressions and the PO sum converges - which it usually does not [13] for non-integrable 
systems in the limit T (or 7) — > 0. Thus, one can emphasize the gross-shell structure in the level 
density using a smoothing width which is much larger than the s.p. level spacing but smaller than 
the main shell spacing (the distance between gross shells) near the Fermi surface. Alternatively, 
finer shell structures can be studied using essentially smaller smoothing widths. 

We note that the POT can also be applied to open systems above the continuum threshold, 
where the level density is dominated by resonances; an example has been given in [53]. As a 
first step towards dynamics, the oscillating parts of nuclear moments of inertia have been studied 
semiclassically in terms of POs in [54, 55]. 

Finally, we mention applications of the POT to finite fermion systems from other domains 
of physics such as mesoscopic semiconductor structures (e.g., quantum dots), metal clusters, and 
trapped fermionic atoms, described in [4] (Chapter 8) and [56, 57, 58, 59, 60]. 



3.2 SEMICLASSICAL CALCULATION OF A NUCLEAR FISSION 
BARRIER 



One prominent feature in the fission of actinide nuclei (isotopes of U, Pu, etc.) is that their fragment 
distributions are asymmetric with a most probable ratio of fragment masses of ~ 1.3 - 1.5 (cf. [12]). 
This is an effect that cannot be described within the LDM which always favours the highest possible 
symmetries. It was one of the big successes of the SCM to explain the mass asymmetry of the fission 
fragments. The fragment distribution is, of course, a result of nuclear dynamics. However, already 
in static calculations of fission barriers, the onset of the mass (or left-right) asymmetry at the 
outer fission barrier was found in SCM calculations with realistic nuclear shell models [61, 62]. On 
the l.h.s. of Fig. 3, we show a schematic picture of the deformation energy of a typical actinide 
nucleus, plotted versus a suitably chosen deformation parameter (see below for a specific choice of 
deformations). The heavy dashed line is the average deformation energy obtained in the LDM; the 
thin lines are the results obtained when the shell-correction energy SE is included. They exhibit 
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the characteristic deformation effects of the shell structure in these nuclei: a deformed ground state 
and the characteristic double-humped fission barrier, split by a second minimum corresponding to 
the fission isomer. The solid line is obtained when only left-right symmetric deformations are used; 
the dashed thin line is obtained when one allows for left-right asymmetric shapes. As we see, the 
asymmetric shapes lower the outer fission barrier considerably. All shapes here are taken to be 
axially symmetric. 

The mass asymmetry in nuclear fission was therefore understood as a quantum shell effect. In 
a detailed microscopical study [63] of the Lund group using the Nilsson model, it was shown that 
those single-particle states which are most sensitive to the left-right asymmetric deformations are 
pairs of states with opposite parity, having the nodes and extrema of their wave functions on parallel 
planes perpendicular to the symmetry axis at and near the waist-line of the fissioning nucleus, as 
shown on the r.h.s. of Fig. 3. Under the effect of the neck constriction one of these s.p. levels, which 
for actinides is just lying below the Fermi energy, is further lowered when the mass asymmetry is 
turned on. As a consequence, the asymmetry leads to a lowering of the total shell-correction energy 
and hence of the outer fission barrier, the LDM part of the energy being much less sensitive to the 
mass asymmetry. (As a historical note, we should emphasize that the essential role of such pairs of 
s.p. states for the mass asymmetry had already been recognized in 1962 by S.A.E. Johansson [30]. 
Unfortunately, at that time the SCM was not available and hence no realistic fission barriers could 
be calculated, as reported in Sect. 2.) 

In this section, we want to show that the POT is able to reproduce this quantum shell effect, 
at least qualitatively, in s semiclassical description using the POT. We will focus here only on the 
gross-shell structure, like that seen in the qualitative picture of the fission barrier in Fig. 3. 

The spheroidal cavity model used in [5] and discussed in Sect. 3.3 does not allow one to describe 
nuclear fission, since an ellipsoidal deformation is not sufficient to yield a finite barrier towards 
fission. In [26, 64, 65], a simple but more realistic "fission cavity model" was used. It consists 
of a cavity with the (c, h, a) shape parametrization that was used both for the LDM and for the 
deformed Woods-Saxon type shell-model potentials in the SCM calculations of [12]. These axially 
symmetric shapes are shown in Fig. 4. The parameter c describes the elongation of the nucleus (in 
units of the radius Ro of a sphere containing N = p 4:7rRl/3 particles, where po is the bulk particle 
density), h is a necking parameter, and q^O describes left-right asymmetric shapes shown by the 
dotted lines. The sequence of shapes with h — a — reproduces the optimized shapes of the LDM 
[33] (see [12] for details). Like in [5], spin-orbit and pairing interactions were neglected in [26, 64, 65] 
and, for simplicitly, only one kind of nucleons (without Coulomb interaction) was used. The only 



parameter in the fission cavity model, the Fermi wave number hp = 12.1/_R , wa s adjusted to yield 
the second minimum at the deformation h = a = 0, c = 1.42 which is that of the fission isomer 
obtained in [12] for the nucleus 240 Pu. This corresponds here to a particle number N ~ 180, i.e., to 
jyi/3 ^ 5 when the spin-orbit interaction is neglected. 

This procedure is justified by the observation that, to a first approximation, the spin-orbit 
and Coulomb interactions essentially lead to a shift of the magic numbers, preserving the relative 
shell structures in the energy shell-correction. This shift can be simulated by a shift of the Fermi 
energy as in [5]. The procedure works, however, only locally in a limited region of deformations 
and particle numbers. The results shown in Fig. 10 below suggest that it is successful in the 
region 1.3 < c < 1.65; the ground-state deformations would, e.g., not be reproduced correctly with 
the same Fermi energy. [Note that, in principle, spin-orbit effects can be included in the POT 
[66, 67, 68] (see also [69] for a short overview). However, in non-integrable systems one is met with 
lots of bifurcations under the variation of the spin-orbit strength [70], which makes the POT with 
spin-orbit interactions very cumbersome. Similarly, the pairing interactions can also be included in 
the POT [71, 72], but this has not been done for nuclear deformations energies so far.] 

In [26], the shortest POs in the (c,h,a) cavity were found to dominate the gross-shell features 
of the double- humped fission barrier. For the deformations around the barriers (c > 1.3), the 
shortest POs are the primitive diagonal and regular polygonal orbits in planes perpendicular to the 
nuclear symmetry axis, situated at the extrema of the cavity shape function (see the l.h.s. of Fig. 
8 below). At the onset of the neck (c = 1.49 for h — a — 0), the orbits in the central equatorial 
plane become unstable with respect to small perturbations perpendicular to the equatorial plane 
and give birth to new stable orbits lying in planes parallel to the equatorial plane. In the restricted 
deformation space with a = 0, these bifurcations are of pitchfork type; they are isochronous due to 
the reflection symmetry with respect to the equatorial plane (cf. [73], Sect. II. B. 2, for a discussion 
of this generic bifurcation type). When the asymmetry a ^ is turned on in the presence of a 
neck, the bifurcation is of a more complicated type. These bifurcations were treated in a uniform 
approximation employing a suitable normal form for the action function (see [26] for details). (Note 
that with respect to small perturbations within the equatorial plane, all these orbits are marginally 
stable, forming degenerate families with degeneracy JC = 1 due to the axial symmetry of the cavity.) 

Before summarizing the results of [26], let us study the general trends of the shell effects obtained 
in the fisson cavity model and try to understand them in terms of the leading POs. 

In Fig. 5 we show a contour plot of the quantum-mechanical shell-correction energy SE calculated 
from the s.p. energy spectrum of the fission cavity model, shown versus the cube-root of the particle 



number iV 1 / 3 and the elongation parameter c along h — a — (white: positive values, gray to black: 
negative values). The horizontal dotted line for N ~ 180 (i.e., k F = 12.1/R ) corresponds to the 
situation where the isomer minimum lies at c ~ 1.4 and the outer (symmetric) barrier is peaked 
around c ~ 1.55, as seen in Figs. 7-10 below. The heavy lines give the loci of constant actions of the 
leading POs (3,l,l)s: meridional triangles; (2,1)EQ: equatorial diameters; (2,1)AQ: diameter orbits 
in planes parallel to the equator plane. We see that these lines follow the valleys of minimal shell- 
correction energy. For the valleys corresponding to the ground-state deformations, the situation is 
like in Figs. 2 (Sect. 1) and 23 (Sect. 3.3 below) obtained for the spheroidal models, but here for 
the more realistic fission cavity model; in all cases the meridional orbits dominate the ground-state 
valleys. The valleys corresponding to the fission isomers, starting around c ~ 1.3, are determined 
by the shortest POs in the planes perpendicular to the symmetry axis: up to c ~ 1.5, these are 
the equatorial orbits EQ; after their bifurcation at c = 1.49, the valleys are seen to curve down 
towards smaller values of iV 1 / 3 , following the constant-action lines of the stable POs in the planes 
parallel to the equator plane (dashed lines, AQ) which have larger semiclassical amplitudes than 
the equatorial orbits (EQ) that for c > 1.49 have become unstable. The fact that the quantum- 
mechanically obtained stability valleys follow the (dashed) lines AQ after their branching from the 
lines EQ is a remarkable quantum signature of the classical bifurcation effect. 

The most striking feature of the gross-shell structure seen in Fig. 5, namely the opposite slopes 
of the ground-state valleys (1 < c < 1.25) and the isomer valleys (c > 1.3), are thus understood 
semiclassically in terms of the opposite deformation dependence of the dominating meridional POs 
in the former valleys and the POs in planes perpendicular to the symmetry axis in the latter valleys, 
respectively. 

These results can be further elucidated by looking at the Fourier spectra in Fig. 6 for the 
five values (from top to bottom) c = 1.1, 1.2, 1.4, 1.5, and 1.6 (all for symmetric shapes with 
h — a — 0). The short arrows underneath the Fourier peaks indicate the lengths of the equatorial 
orbits: diameter (2,1)EQ and its second repetition 2(2,1)EQ, triangle (3,1)EQ, etc., and (for c = 1.6) 
the corresponding orbits AQ in the planes parallel to the equator plane. The long arrows correspond 
to the meridional orbits: triangle (3,l,l)s and quadrangle (4,1, l)s. We see that for the small 
deformations c = 1.1 and 1.2, the meridional orbits have the strongest amplitudes and hence 
dominate the shell structure in yielding the ground-state deformation valleys seen in Fig. 5. The 
equatorial orbits EQ and their bifurcated partners AQ have the largest amplitudes for c = 1.4 — 1.6, 
which explains their dominance in yielding the isomer valleys. 

Since the trace formula (12) for the shell-correction energy SE contains the squares (h/t po ) 2 of 



the inverse periods of the POs, apart from the amplitudes A po of the semiclassical level density (11), 
there is a natural suppression of longer orbits contributing to SE. This ensures the convergence of 
the PO sum, particularly in non-integrable systems (like the one considered here) where the PO sum 
for the level density (11) usually does not converge [13]. This suppression is particularly effective 
amongst orbits with comparable amplitudes A po . It explains why already at c = 1.4, where the 
meridional orbits (3,l,l)s and (4,1, l)s still have similar amplitudes as the EQ orbits, 2 the latter 
dominate the shell structure (by a factor ~ 4 in case of the EQ2 orbit), as suggested by Fig. 5. 

In Fig. 6 we have marked some of the peaks around 7 < L < 8.5 for c = 1.4 and around 
6.5 < L < 8 for c = 1.5 and 1.6. They correspond to orbits born from the equatorial orbits in 
period-doubling bifurcations; some of them are 3-dimensional orbits. Similarly, there are many other 
peaks at L > 8, some of which correspond to orbits born in higher m-tupling (m > 3) bifurcations. 
The contributions of all these orbits to the gross-shell structure is, however, practically negligible 
due to their long periods. They have therefore not been included in the results presented below. 
They might, however, become noticeable in POT calculations with higher resolution of the shell 
structure (see the discussions in Sects. 3.3 and 4). 

We should also recall the fact that in realistic SCM calculations, the pairing interactions are 
known to reduce the amplitude of 5E by up to ~ 30% (see, e.g., [12]). In the POT, the pairing 
effects [71, 72] yield, indeed, an extra smoothing factor in the semiclassical amplitudes, which further 
suppresses the contributions of longer orbits. 

Fig. 7 shows a comparison of the level density in the fission cavity model, taken at the Fermi 

wave number kp = 12.1/ Rq and plotted versus c along h = a = in the region covering the isomer 

minimum and the outer fission barrier. [As noted above, the model without spin-orbit interaction 

and one fixed Fermi energy does not work at all deformations. For this reason the deformations 

c < 1.25 are not shown.] The solid line is the quantum- mechanical result and the dashed line 

the result of the semiclassical trace formula (with uniform approximation [26] for the bifurcation 

at c = 1.49) including all primitive POs with up to 20 reflections at the boundary in the planes 

perpendicular to the symmetry axis. Both results have been coarse-grained by a Gaussian averaging 

over the wave number k with a width 7 = 0.6/i?o that emphasizes the gross-shell structure. We 

see that the agreement is nearly perfect. This figure also demonstrates that already a very simple 

cavity model with one degree of freedom (c, keeping h = a = fixed) is capable of yielding the 

2 The strong Fourier peak near L/R ~ 6.3 for c = 1.4 in Fig. 6 contains the combined amplitudes of the meridional 
quadrangle (4,l,l)s and the second repetition of the equatorial diameter orbit, 2(2,1)EQ. Although the two cannot 
be disentangled, we estimate that both these orbits have comparable amplitudes. 



main features of the double-humped fission barrier. 

Let us now look at the influence of left-right asymmetric shapes with a ^ Oon the shell-correction 
energy and, in particular, on the height of the second fission barrier. In Fig. 8, the semiclassical 
result of SE is shown in a perspective view as a function of elongation c and left-right asymmetry 
a, taken along h = in the region of the isomer minimum and the outer fission barrier. We see how 
the outer fission barrier is lowered for left-right asymmetric shapes. Instead of the heigher barrier 
obtained for left-right symmetric shapes with a = (arrow labeled "symm."), the nucleus can go 
towards fission over a lower saddle when asymmetric shapes are allowed (arrow labeled "asymm."). 
To the left, we see the shapes corresponding to the three points A (fission isomer), B and C (along 
the asymmetric fission path) in the deformation energy surface. The vertical lines indicate the 
planes in which the POs are situated (solid lines for stable and dashed line for unstable POs). 

The instability of the outer fission barrier towards left-right asymmetric deformations, known 
from the quantum-mechanical SCM calculations [61, 62], can thus be described semiclassically using 
the POT, indeed. Hereby only the shortest primitive POs are relevant due to the fast convergence 
of the PO sum for the semiclassical shell-correction energy SE, as discussed above. 

This convergence is also demonstrated in Fig. 9. It shows the Fourier transform of the level 
density obtained for the asymmetric shape with c=1.5, h=0, a=0.12, corresponding to the point B 
shown in Figs. 8 and 10 below, coarse-grained as in Fig. 7. The insert depicts the corresponding 
shape. The solid line is the quantum- mechanical and the dashed line the semiclassical result; the 
two are seen to agree very well. We notice that the most significant peaks correspond to the lengths 
L/R of the primitive orbits (p, 1)AQ with p — 2, 3, 4, and 5 reflections. Some of the small peaks 
at L > 6R correspond to their second repetitions, the other come from meridional orbits. These 
do, however, not affect the gross-shell structure (cf. also Fig. 4 in [64]). 

In Fig. 10, we compare the old quantum-mechanical results of SCM calculations [12] with realistic 
deformed Woods-Saxon potentials [11] with the semiclassical POT results using the present simple 
fission cavity model. Shown are contour plots of SE versus c and a for two values of the neck 
parameter h. We see that the semiclassical results (r.h.s.) reproduce the gross-shell structure of 
the quantum results (l.h.s.) very well. The correct topology is obtained, displaying the lowering 
of the outer barrier for left-right asymmetric shapes. Also, the amplitudes of the shell effects on 
both sides are comparable, which justifies our calculating only the gross-shell structure using the 
shortest periods on the semiclassical level. Of course, a detailed quantitative agreement cannot be 
expected for the two calculations using such different potentials as the sophisticated smooth Woods- 
Saxon potential including pairing, spin-orbit, and Coulomb interactions on one side (left), and the 



simple fission cavity model without these extra interactions on the other side (right). The more 
gratifying is the good overall good qualitative agreement of the gross-shell structure. This agreement 
demonstrates, by the way an experience made from quantum-mechanical SCM calculations using 
realistic nuclear shell model potentials: the gross-shell features of the fission barriers are much less 
sensitive to the radial dependence of the potentials than to its deformation. Hence the success of a 
simple cavity model that is very schematic, but uses the realistic c, h, a deformations. 

The white dashed lines in the r.h.s. panels of Fig. 10 show the loci of constant classical actions 
S'po of the leading POs. They follow exactly the valleys of minimal energy in the (c, a) planes 
which define the adiabatic fission paths. Thus, as it was already observed in [5] and seen in Fig. 
5, the condition for minimizing the shell-correction energy is semiclassically given by a least-action 
principle: 5Spo = 0. 

We should emphasize that in Figs. 8 and 10 only the shell-correction energy 5E is shown. The 
complete fission barrier is obtained by adding its smooth LDM part, according to (1), which for 
240 Pu in the (c, h, a) parametrization occurs [12] at c ~ 1.45 — 1.5, h — a — 0. Since the LDM 
barrier is rather smooth around its maximum, the relative heights of the isomer minimum and the 
outer barrier are not affected much by it. However, for c > 1.6 the LDM barrier is already going 
steeply down. Therefore in the total energy, the minimum around c ~ 1.65 — 1.7 seen in Figs. 8 and 
10 along h = = a = (like also in Fig. 7) vanishes in the steep slope of the total fission barrier, 
as shown schematically in Fig. 3. 

It is also interesting to note, as demonstrated explicitly in [65], that the quantum-mechanical 
probability maxima of those s.p. states which microscopically are responsible [63] for the asymmetry 
effect in the SCM approach (see the schematic plot on the r.h.s. of Fig. 3) lie exactly in the planes 
perpendicular to the symmetry axis that contain the classical POs. This constitutes a nice quantum- 
to-classical relationship. It was also shown in [65] that the classical dynamics of the nucleons with 
small angular momenta L z is more than 90% chaotic in the region of the outer barrier. A very small 
phase-space region of regular motion is thus sufficient to create the shell effect that leads towards 
the asymmetric fission of the nucleus! 

We emphasize once more that the fission cavity model, in its present form without spin-orbit 
and Coulomb interactions, is not suitable for predicting fission barriers for a larger range of nuclear 
isotopes and deformations. The present semiclassical calculation should be taken as a model study 
of a typical actinide fission barrier, demonstrating that the POT in principle is capable of explaining 
the existence of a double-humped barrier, and also the onset of mass asymmetry around the outer 
barrier, in terms of a few short classical periodic orbits. It was in no way meant as a substitute 



for the quantum-mechanical SCM for calculations of static fission barriers. Its aim was rather to 
provide, as suggested by the late Strutinsky, a qualitative physical understanding of a sophisticated 
quantum shell effect by means of simple classical pictures. 



3.3 SPHEROIDAL CAVITY 



As one of the main subjects in this report, we apply the semiclassical theory to a 3D spheroidal 
(prolate) cavity [5, 19, 20], which may be taken as a simple model for a large deformed nucleus or a 
(highly idealized) deformed metallic cluster [4]. We shall investigate the role of orbit bifurcations in 
the shell structure in connection with the superdeformation (corresponding to the fission isomer of 
heavy nuclei). Although the spheroidal cavity is integrable, it exhibits all the difficulties mentioned 
above (i.e., bifurcations and symmetry breaking) and therefore gives rise to an exemplary case study 
of a nontrivial 3D system. We apply the ISPM (see sect. 3.1) for the bifurcating orbits and succeed 
in reproducing the superdeformed shell structure by the POT, hereby observing a considerable 
enhancement of the oscillation amplitudes of the level density and the energy shell corrections near 
the bifurcation points and for larger deformations. 

3.3.1 CLASSICAL DYNAMICS 

We consider a cavity system having a spheroidal shape boundary with the lengths of minor and 
major axis a and b, respectively, with a conserved volume a 2 b = R 3 . The spheroidal deformation is 
therefore described by a single parameter, which we take to be the axis ratio rj — b/a. The classical 
dynamics of particles in the spheroidal cavity can be solved in terms of action-angle variables related 
to the elliptic coordinates k — v,u, ip in this Hamiltonian system. The Hamiltonian, H(I V ,I U ,I V ), 
depends only on the partial actions [5, 20], 



where I\ and I2 are new single-valued separation constants, introduced for convenience instead of 
the two remaining free actions for the fixed energy by using the energy conservation of the particle 
motion: E = H(I U , I v , I^,), and the boundaries of the spatial part of the phase-space occupied 
by classical trajectories are shown explicitly in the integration limits. (These two constants are 
identical to o\ and a 2 in the notation of Refs. [5, 20]). The periodic-orbit (SPM) equations (7) for 
the stationary points I{ and I^ can be written as a resonance condition of commensurability of the 




(18) 



partial frequencies: 

uj v : u u : Up = n v : n u : , with cu K = dH/dI K . (19) 

For definiteness, we shall consider only prolate spheroidal deformations, rj > 1, following the sug- 
gestions in [5]. 

The POs in the spheroidal cavity can be classified as follows: a) meridional (elliptic, 2DE) 
orbits, denoted by M (n v ,n u ) (n v > 2n u ), and b) equatorial (EQ) orbits, denoted by M (n v ,n v ) 
(n v > 2n v ). Hereby M is the repetition number of the primitive orbits which are characterized by 
the positive integers n v = number of corners and n u ,n^ = winding numbers), see Figs. 11 and 12, 
respectively. Both these types of orbits exist at all deformations rj > 1. c) There are also specific 
families of orbits which exist only at specific deformations larger the critical (bifurcation) value, 
i] > rjbif, where they are born from the EQ POs in pitchfork-like bifurcations. Examples of such 
bifurcated orbits M(n v , n u , n v ) are the meridional hyperbolic orbits (2DH), e.g., the butterfly (4,2,1) 
(see the fourth orbit in Fig. 11), which bifurcates from the doubly repeated EQ diameter 2(2,1) (cf. 
the first orbit in Fig. 12) at the deformation r] bi f = V% and exists for all deformations rj > \[2 along 
with the diameter 2(2,1). Another important example is the bifurcation of the equatorial star orbit 
(fourth orbit in Fig. 12) at rj = 1.618... where the new 3D orbits (5,2,1) emerge from the parent 
orbit (see the first orbit in Fig. 13); again both exist at all larger deformations. One also finds a 
bifurcation of the doubly repeated EQ triangle 2(3,1) which bifurcates at 77 = 1.732, whereby the 
new 3D orbits (6,2,1) emerge. 

For the following it is useful to note the classical degeneracies /C of the above-mentioned POs 
[15, 18, 20, 21]. The meridional (planar) orbits (denoted by "2D" henceforth for both elliptic 2DE 
and hyperbolic 2DH POs) and all 3D orbits have the degeneracy /C = 2, so that the largest phase- 
space volume is covered by these PO families. The POs in the equatorial plane (denoted by "EQ") 
have the degeneracy K = 1 (except at the critical bifurcation points where they are degenerate with 
the newborn orbit families having /C = 2). There is also an isolated (JC = 0) diametric orbit along 
the symmetry axis for rj > 1, but its contribution can be neglected for sufficiently large deformations 
as compared to all other orbits. 

In the following subsections, we will show more explicitly that bifurcations of POs are not just 
a technical obstacle which appears in the semiclassical POT, but that they leave signatures in the 
quantum-mechanical results of the system and thus cause physical phenomena like in Sect. 3.2. 



3.3.2 TRACE FORMULAE AND SHELL STRUCTURE ENHANCEMENT 

The contributions of the above periodic orbits are in the following separately considered according 
to their degeneracies K discussed above. According to Eq. (11), the oscillating part of the level 
density is thus given semiclassically by 

5g scl (E) = 5g { £{E) + 8g™{E) + 8g^(E), (20) 
5g^\E) = ReJ>£) exp (ikL po - i\ofJ - i<f>f } ) . (21) 

po 

For the explicit analytical expressions of the amplitudes , constant phases and (f>^\ see 
Ref. [20]. 

Figs. 14 and 15 show the characteristic enhancement of the semiclassical amplitudes near the 
bifurcations of the EQ orbits (5,2) and 2(3,1), respectively. As clearly seen from these Figures, the 
ISPM amplitude (solid lines) Aeq has a pronounced but finite maximum instead of the divergence 
of the SSPM amplitude (dashed lines) at the corresponding critical points 1.618... and 1.732... . 
The SSPM amplitudes A^d for the bifurcated 3D families (5,2,1) and (6,2,1) have discontinuities 
at these deformations. Their ISPM amplitudes are continuous curves with an enhancement in the 
form of a rather wide maximum slightly to the right of the critical deformations. The relative 
enhancement near the bifurcations is of order fT 1 ^ 2 [or (k F L po ) l l 2 at the Fermi energy E = E F \. As 
noted in sect. 3.1, it is due to one more exact integration in the improved SPM, as compared to the 
standard SPM. This enhancement occurs around the critical deformation (on both sides) for the 
parent orbits and to the right of the bifurcation point for the newborn POs. For any integrable and 
non-integrable systems, one has an increase of the classical degeneracy of the parent orbits locally 
at the bifurcation point. The newborn POs form families having the larger degeneracy (/C = 2) at 
the critical point and all larger deformations, as noted above. This discrete behavior of the classical 
degeneracies is "smeared out" semiclassically in the ISPM over the regions around the bifurcation 
point (or slightly on its right) where the enhancement is found in the amplitudes of Figs. 14 and 15 
for the general reasons discussed in Sect. 3.1. 

We observe also some smaller oscillations in these amplitudes further away from the bifurcation 
points. They can be neglected with respect to the enhancement at the order of the semiclassical 
parameter discussed above. In any case, to improve the SPM with respect to its standard version 
(SSPM), we have to use finite limits in the integration over a phase-space volume occupied by CT 
trajectories as given in the more exact expression (3) for the level density (before applying further 
approximations to the trace of the semiclassical Green function). 

Within the simplest ISPM, using the expansion of the action phase near the stationary point up 



to second-order terms like in (9) ("2nd-order ISPM"), one can extend the integration limits to ±00 
for deformations asymptotically far from the bifurcations in order to obtain the SSPM asymptotics. 
For more exact higher-order expansions of the action (8), for example up to third-order terms ("3rd- 
order ISPM"), the convergence for the extension of the integration limits to an infinite region can 
be sufficiently good, even near the bifurcations, see [20]. As shown from the comparison of the 2nd 
and 3rd order ISPM approaches in Figs. 11 and 12 in [20], the convergence of this expansion around 
the bifurcation points is rather fast and, therefore, it can be stopped at the 2nd order ISPM, but 
with finite limits. 

3.3.3 COMPARISON WITH QUANTUM RESULTS 

We now compare numerically our semiclassical results for the spheroidal cavity with exact quantum- 
mechanical results. In Figs. 16 and 17 we show the oscillating part of the level density Sg 1 (k) versus 
the wave number k for six charecteristic values of the deformation parameter r]. Both the quantum- 
mechanical (solid lines) and the semiclassical (dotted lines) level densities have been averaged over a 
Gaussian width of 7 = 0.3/i?, in order to emphasize some finer shell structure besides the gross-shell 
oscillations. In Figs. 18 and 19 we show the results for the energy shell-correction 8E. The quantum 
results (solid lines) are obtained by the shell-correction method, see (2), from the quantum spectrum 
of the spheroidal cavity. The semiclassical results are obtained using (12) which, in analogy to (20), 
becomes 

SE scl = Re £ (f) A$\E F ) exp (ik F L po - ?-ofJ - v^f) , (22) 

K,po V W V Z ' 

where the sum is over the same three orbit families as in (20) and t po = mL po /p F are their periods, 
taken at the Fermi momentum p F = hk F = \/2mE F which is fixed for each particle number N 
through (13). 

A close examination of contributions of individual POs at the deformation r] — 1.2, see Figs. 
16 and 18, shows that a good agreement between the semiclassical (ISPM or SSPM) and quantum 
calculations is obtained by including only short (period-one) elliptic 2D (2DE) and EQ orbits. More 
specifically, we find a rather good convergence of the PO sums by using the 2DE orbits with n v < 12, 
n u — 1 and the EQ orbits with maximum vertex and winding numbers Mn v = 14 and Mn^ = 1, 
respectively. For small and medium deformations up to rj ~ 1.4, one has a clear dominance of the 
gross-shell structure by the shortest meridian and equatorial POs, similar to the situation in the 
fission cavity model (see Fig. 5) up to c ~ 1.4 (see also the discussion below in connection with 
Figs. 21 and 23). 

At 77 = \/2, there occurs a bifurcation of the doubly repeated EQ diameter 2(2,1) with the emer- 



gence of the butterfly orbits (4,2,1)2DH, and the improved SPM becomes important in comparison 
with the quantum results, see the second panels of these figures. However, even at these deforma- 
tions (77 < 1.4), as shown in [20], these longer orbits do not influence much the gross-shell structure 
due to the leading period-one orbits mentioned above. Noticeable ISPM contributions of the bifur- 
cating (4,2,1)2DH orbits are found in the lowest panels of Figs. 16 and 18 for r\ = 1.5, in agreement 
with the enhancement of their amplitudes to the right of their bifurcation point. At r\ = 1.5, these 
contributions become comparable in magnitude with respect to those of the period-one orbits, al- 
though they have similar periods, and they become even larger with increasing deformation and 
particle number. Thus, for 77 > 1.5 one has a transition from period-one to period-two dominating 
PO contributions; the latter are strongest at i] = 2 and beyond. 

Figs. 17 and 19 show the cases of superdeformed states exactly at the bifurcation points rj = 
1.618..., 1.732... and 2.0. The quantum results are here nicely reproduced by the semiclassical ISPM 
results (like for r\ = 1.414... above) including the dominant contributions from the bifurcating orbits. 
Their important role is illustrated in Fig. 20 for the same superdeformations rj > 1.6, where only 
a few bifurcated period-two orbits are used in the ISPM results, which are again compared to 
the quantum results. It is seen from this Figure that the main features (periods and amplitudes) 
of the shell and supershell structures are in fairly good agreement which becomes the better the 
larger particle number. Some of the remaining discrepancies are still due to the (much smaller) 
contributions of period-one orbits and some other missing period-two orbits. 

The importance of the bifurcating period-two orbits at these larger deformations is further 
stressed in the following two figures. In Fig. 21, the thin solid and dotted lines are again the 
quantum and semiclassical ISPM results, like in Figs. 18, 19 and 20, showing a good agreement. 
The heavy line shows the ISMP result which is obtained if one uses only the shortest (period-one) 
orbits. We see that with increasing values of r] and iV 1//3 , the resulting gross-shell structure does 
not reproduce the basic features mentioned above. At the super deformation i] = 2, the shortest 
alone orbits fail, in fact, completely to reproduce the gross-shell structure. At rj = 1.732..., the shell 
structure is dominated mainly by the (4,2,1)2DH and (5,2,1)3D orbits, which have emerged from 
bifurcations at rj^f = 1.414... and 1.618..., respectively, and by the 2(3,1)EQ and (6,2,1)3D orbits 
bifurcating exactly at rjuf = 1.732.... At 77^/ = 2, the (6,2,1)2DH orbits bifurcate from the 3(2,1)EQ 
orbits. Besides these and the above-mentioned POs, also other new bifurcating period-three orbits 
like (7,3,1) and (8,3,1) 3D - besides the 3(2,1)EQ orbits which have enhanced amplitudes at rj — 2 
as usually - are determining the semiclassical shell-correction energies. 

These findings are also transparently supported in Fig. 22 by the quantum Fourier spectra at the 



same deformations as above (see a more detailed discussion of the Fourier spectra below). Many 
high peaks in the superdeformed region 77 = 1.6 — 2.0 remarkably confirm that the period-two 
POs are mainly responsible for the shell structure here and dominating above the period-one POs, 
especially clearly so at 77 ~ 1.7 — 2.0. 

We see, therefore, that in the superdeformed region, the period-one orbits alone cannot explain 
the shell structure. The reason why the bifurcated period-two orbits become so important lies in 
their enhanced semiclassical amplitudes (cf. Figs. 14 and 15 above) which significantly overcome 
the suppression due to the l/£~ factor in (22). One important ingredient of this enhancement is 
the larger phase-space volume covered by the newborn orbits due to their degeneracy JC = 2 (see 
Sect. 3.1 for a more general discussion of the origins of the enhancement). 

In Figs. 23 and 24 we plot the quantum-mechanical shell-correction energy SE as function of 
deformation 77 and particle number iV 1//3 in order to study its shell structure. The basic gross- 
shell structure of the energy valleys of minima is nicely described by the constant-action curves 
of the leading shortest POs at deformations smaller than the superdeformed shapes, as shown 
in Fig. 23. The main features are rather similar to those seen in Fig. 2 for the Woods-Saxon 
potential with spheroidal shapes, taken from [5], and also similar to Fig. 5 for relatively small 
deformations in the more realistic fission cavity model. In the region of small deformations, the 
dominant period-one orbits are the meridional 2DE orbits (3,1,1) and (4,1,1) whose constant-action 
lines follow precisely the down-going ground-state valleys in all three models. At larger intermediate 
deformations 1.3 < 77 < 1.5, the slopes of the upgoing second-minimum valleys are roughly followed 
by the constant-action lines of the period-one orbits (2,1)EQ and (3,1)EQ, which is similar again 
as in Fig. 5 for the fission cavity model for 1.25 < c < 1.5. However, for 77 > 1.6 the secondary 
minimum valleys have a trend to bend down towards smaller values of N 1 ^ 3 with increasing 77 , 3 
which is not explained by the primitive EQ orbits. This trend can only be explained by the 
influence of the longer orbits which bifurcate from the higher repetitions of the EQ orbits. As 
seen from Fig. 24, the constant-action lines of the newborn orbits (4,2,1), (5,2,1) and (6,2,1) do, 
in fact, follow the downbending of the second-minimum valleys very nicely, especially for 77 > 1.7 
in good agreement with Fourier peaks in Fig. 22. In the region where the downbending begins, 
77 = 1.4 - 1.5, the dominance of the 2(2,1)EQ and the newborn 2DH orbits born at 77^/ = 
reflects the enhancement of their semiclassical amplitudes at and beyond the critical deformation. 
This confirms the conclusions drawn above in connection with Figs. 21 and 22. In the region of 
superdeformations 77 > 1.6, the newborn period-two and period-three orbits dominate the gross- 
3 Note that this trend can also be seen in Fig. 2 for the diffuse Woods-Saxon potential with spheroidal deformations. 



shell structure of the shortest (period-one) EQ orbits completely, as clearly seen also from Fig. 
22. 

[Note that similar properties of the superdeformed shell structure owing to the amplitude en- 
hancement of the dominating bifurcated period-two (and higher period) orbits were shown in more 
detail, including the comparison of the maps of the energy shell corrections 5E(rj,N) for quantum 
and semiclassical calculations (with and without bifurcating POs), in Fig. 20 of [18]. In particular, 
the ISPM results for SE were in nice agreement with the quantum results at all critical characteristic 
deformations. The bifurcating period-two orbits were important to reproduce the superdeformed 
shell structure at rj > 1.6 like for the present spheroid cavity] 

We point out that the type of the bifurcations causing the downbending of the isomer valleys in 
the fission cavity model, seen in Fig. 5, is a different one: there it is the period-one bifurcation of 
the primitive EQ orbits at c = 1.49 which gives birth to the AQ orbits (also of period one) whose 
constant-action lines follow correctly the downbending for c > 1.5. It is a subject of our ongoing 
research to investigate how bifurcations of further period-one and period-two orbits in the fission 
cavity model will affect finer details of the shell structure and, in particular, the "quantization" of 
the various local minima along the isomer valleys (see also the discussion in Sect. 4). 

Note that in the axially-symmetric deformed 3D harmonic oscillator, discussed in [5], a smooth 
upgoing behaviour of the valley of minimal SE is more pronounced due to the shortest EQ orbits 
along their constant-action lines (following the gross-shell structure). It is accompanied by deeper 
local minima along the valleys (see Fig. 2d in [5]). They are associated with the amplitude en- 
hancement owing to the longer-3D orbit bifurcations causing finer shell structure details. These 
local minima are related to the bifurcations of the simplest 3D orbits, for instance (3,2) at the de- 
formation i] = u±/u z = 1.5 and (2,1) at 2.0, where oj±_ and u z are the HO frequencies perpendicular 
and parallel to the symmetry axis. The bifurcation amplitude enhancement near the intermediate 
deformation rj ~ 1.5 there does not destroy the leading gross-shell structure valleys, like in the two 
cavity models discussed here. However, the decreasing slopes at superdeformations r\ 2 are much 
less pronounced in the deformed HO model. We emphasize the important difference of the cavity 
models with respect to the deformed HO that the enhancement of the amplitudes of the bifurcating 
orbits occurs for all deformations larger than the bifurcation point, see [5]. We may also point out 
that the downgoing ground-state deformation valleys, due to the meridional orbits in the present 
cavity models, are absent in the harmonic oscillator and, in fact, degenerated into one spherical 
point r] = 1. Each of the bifurcating 3D orbits appears and exists only at one deformation in 
the HO, whereas they exist in both cavity models at all deformations larger than their bifurcation 



points. 

The above discussed properties of the gross and fine shell structure related to the shortest and 
the bifurcating period- two orbits in the spheroidal cavity are nicely confirmed also by Fig. 25, where 
we show the modulus of the Fourier transform \F(L)\, see Eq. (14), of the quantum spectrum of 
the spheroidal cavity, plotted versus the orbit-length variable L and deformation 77. It exhibits 
clear quantum signatures of the amplitude enhancement owing to the newborn orbits above their 
bifurcation points (marked by black dots). Besides the enhancement of \F(L)\ at the PO lengths 
(3,1), (4,1) (and others) at the spherical shape (rj = 1), we observe the enhancement also above 
the bifurcation of the butterfly orbit at 77 = \/2 and the superdeformations 1.618..., = 1.732... 
etc., with the specific PO lengths L = £(5,2,1), ^6,2,1), etc. The enhanced Fourier signals, which 
correspond to the highest peaks in Fig. 22, follow precisely the various lines in Fig. 25 which show 
the lengths of the classical POs (as explained in the insert on the upper right). This is again a nice 
quantum-to-classical correspondence, showing that classical periodic orbits and their bifurcations 
do have a physical relevance in the quantum spectrum of this system. 

It is worth noting that the semiclassical origins of the deformed shell structures discussed above 
are related with bridge orbits [56] , at least for potentials with finite diffuseness which can be approx- 
imately described by the power-law potential discussed in the following subsection. In a power-law 
potential with spheroidal shape, we have strong enhancement of shell effects due to the bridge orbit 
bifurcation (n : m) which emerges from a period-n-tupling bifurcation of equatorial orbits and sub- 
merges into a symmetry-axis orbit through a period-m-tupling bifurcation, forming the "bridge" 
between those POs, (see [56] for details). The meridional 2DE orbits (n, 1, 1) responsible for the 
deformed shell structure at small deformations correspond to (1:1) bridges, and the meridional 2DH 
and 3D orbits (n, 2, 1) responsible for the superdeformed shell structures correspond to (2:1) bridges, 
and so on. In the next subsection, we discuss the power-law potential with spherical shape, which 
also shows bifurcation enhancement effects at certain values of the diffuseness parameter. 

3.4 RADIAL POWER-LAW POTENTIALS 

The idea of [28] is that the spherical Woods-Saxon (WS) potential, known as a realistic mean-field 

potential model for spherical nuclei and metallic clusters, is nicely approximated (up to a constant 

shift and without the spin-orbit term) by much a simpler power-law potential 4 which is proportional 

4 Note: The potential power parameter a here should not be confused with the asymmetry deformation parameter 
of Sect. 3.2. 



to a power of the radial coordinate r a 



V(r) = U(r/R) a . (23) 

With a suitable choice of the parameters U and a, the approximate equality 

Vws(r) ~ V ws (0) + U(r/R) a 

holds up to r < R, where R = r^A 1 ^ with r$ = 1.2 fm represents the nuclear radius for a given 
mass number A. Thus one finds nice agreement of the quantum spectra up to Fermi energy. 

The potential (23) includes the limits of the harmonic oscillator (a = 2) and the cavity (a — > 
oo); realistic nuclear potentials with steep but smooth surfaces correspond to values in the range 
2 < a < oo. The advantage of this potential is that it is a homogeneous function of the coordinates, 
so that the classical equations of motion are invariant under the scale transformations: 

r s 1/a r, p -> s 1/2 p, t s 1/2 ~ 1/a t for E sE . (24) 

Therefore, one only has to solve the classical dynamics once at a fixed energy, e.g. E = U; the 
results for all other energies are then simply given by the scale transformations (24) with s = E/U. 
This highly simplifies the POT analysis [28, 56]. 

With this, we are able to apply the POT to more realistic potentials, accounting for a finite 
diffuseness of the potential along the radial variable. Note that the definition (23) can also be 
generalized to include deformations (see, e.g., [56]). 

3.4.1 PERIODIC ORBITS IN THE RADIAL POWER-LAW POTENTIALS 

The definition (23) can be used in arbitrary dimensions, as long as r is the corresponding radial 
variable. In practice, we are interested only in the 2D and 3D cases. The spherical 3D and the 
circular 2D potential models have common sets of periodic orbits. They are labeled by integers 
MlfirjU^), where n r and are mutually commensurable integers representing the number of 
oscillations in radial direction and the number of rotations around the origin, respectively, and M 
is again the repetition number. For the harmonic oscillator (a = 2), all the classical orbits are 
periodic with (degenerate) ellipse shapes. By slightly varying a away from 2, the diameter and 
circular orbits remain as the shortest periodic POs. With increasing a, the circular orbit and its 
repetitions cause successive bifurcations, generating various new periodic orbits (n r ,n^) (n r > 2n v ) 
(see Appendix A2 for details). Fig. 26 shows some of the shortest POs. The only difference between 
the periodic orbits in the 3D and 2D systems is their degeneracies. The shortest PO is the diameter 



which has the degeneracy /C = 1 in the 2D system and /C = 2 in the 3D case at a > 2. Other 
polygon-like orbits have, respectively, /C = 1 and /C = 3 at a > a bi j where a bi j is the bifurcation 
value (A14). The circular orbit, having maximum angular momentum, is isolated (/C = 0) for the 
2D system and has K, = 2 in the 3D case (except at the bifurcation points). 

3.4.2 TRACE FORMULA FOR 2D CIRCULAR POTENTIAL 

In the 2D circular potential, it is useful to take the angular momentum £ as one of the integration 
variables in the phase-space trace formula (3). The integration limits £ = ±£ c , given in (A9), 
correspond to the circular orbits moving anticlockwise (+) and clockwise (— ). 

For the contributions of the one-parametric families of the diameters and polygon-like POs, the 
integration over the cyclic angle variable ip conjugate to £ is easily carried out, and the remaining 
integration over £ is done within the simplest ISPM with finite limits — £ c < £ < £ c . Thus, one 
obtains the ISPM trace formula in terms of error functions, as explained in sect. 3.1. Extending 
the integration limits to ±oo, one will arrive at the asymptotic Berry- Tabor trace formula [16], like 
in the cavity model. It is also worth noting that the limits a — > oo and a — > 2 of the derived ISPM 
trace formula successfully reproduce their asymptotic limits obtained previously for the circular 
billiard (disk) and the 2D isotropic harmonic oscillator, respectively. 

The contribution of the isolated (/C = 0) circular orbits are derived by adopting the ISPM for two 
integrations over suitable (radial phase-space) variables. It yields a product of two error functions 
as shown in (9). Extension of the finite integration limits in both error functions to infinity (i.e., 
to ±oo after their transformation to the corresponding Fresnel integrals of positive or negative real 
arguments) leads to the (SSPM) Gutzwiller amplitudes for the considered isolated POs [13]. 

The total semiclassical level density (5) is thus given by 

9sci(E) = g TF {E) + 6gV(E) + 6g<?J(E) (25) 
with g (E) gxF(E) of the Thomas- Fermi model and 



S9%{E) = ReJ^A^iE) exp ^(E) - |«#> - iff 



po 



(26) 



The analytic expressions for the amplitudes A p ^ and phases and 0^ are given in Appendix 
A3. Using the scale invariance (24), one may factorize the action integral S po (E) as 

Spo(E) = I pdr= (^] I pdr = her po . (27) 



po(E) \U y Jpo(E=U) 



In the last equation, we define dimensionless scaled energy e and scaled period r po by 



£\ 5 



= \l pdr. (28) 
n J vo (E=u) 



To realize the advantage of the scaling invariance, it is helpful to use the scaled energy (period) 
in place of the corresponding original variables. For the harmonic oscillator one has a — 2, and 
the scaled energy and period are proportional to the unsealed quantities. For the cavity potential 
(a — > oo), they are proportional to the momentum p and length L po , respectively. The energy-scaled 
level density is written as 



= d 4^\E) = Re5>£)(e)exp 



de 

po 



i£r po - jv P o ] ~ i<t>f ) 



(29) 



This simple form of the phase function enables us to make use of the Fourier transformation tech- 
nique again. The Fourier transform of the semiclassical energy-scaled level density with respect to 
the scaled period r is given, like in Eq. (14), by 

F{t) = J de g(e)e t£T ~ F (r) + ]T A po 6(r - r po ) , (30) 

po 

which exhibits peaks at periodic orbits r = r po . F (t) represents the Fourier transform of the 
smooth TF level density and has a peak at r = related to the direct (zero-action) trajectory, as 
discussed in Sect. 3.1. Thus, from the Fourier transform of the energy-scaled quantum- mechanical 
level density, 

r i -+— 

F (r) = [ E 5 ( £ ~ £ «) el£Td£ = E e " nT ' £ - = 2 " » ( 31 ) 

one can directly extract information about the classical PO contributions. 

3.4.3 AMPLITUDE ENHANCEMENT AND COMPARISON WITH QUANTUM RESULTS 

Fig. 27 shows the typical enhancement phenomena owing to bifurcations of POs, e.g. for the case 
of the birth of the triangle-like PO (3,1) from the simplest (primitive) circular orbit (1,1). The 
scaled amplitudes |^4 SCiPO | are presented as functions of the power parameter a. In Fig. 27, the 
enhancement of the diameter amplitudes \A 8Cj (2, 1)| and those |^4. SCj j for the primitive circular 
orbit (1, 1) in the HO limit a — > 2 are clearly shown. The second sharp peak in the circular PO 
amplitude A sc , (1,1) near the bifurcation point a = 7 is related to the bifurcation of the PO (3, 1) 
from the parent orbit (1,1). Note that a rather wide maximum appears in the ISPM amplitude 
|Asc(3,i)| of the triangular- like orbit to the right of the bifurcation a = 7 in Fig. 27, even larger than 



for 3D orbits in the spheroidal cavity (see Sect. 3.3). [The third maximum in \A SC , (i,i)| for a birth 
of the simplest quadrangle-like PO (4,1) at the bifurcation a = 14 of the same parent-circular orbit 
is also seen, because the convergence of the ISPM triangle (3,1) amplitude \A SC , (3,1) | to its SSPM 
asymptotics is realized within such a very wide region up to about a = 20.] The ghost part of the 
\A SC , (3,i)l curve for a to the left of the bifurcation point a = 7 is suppressed by rapid oscillations 
of the phase of a complex amplitude A SCt (31) with a frequency going to zero in the HO limit. It is 
similar like for all 3D PO ghosts in the spheroidal cavity in Sect. 3.3 (see also [20]), and for such 
one-parametric family POs in the integrable version of the 2D Henon Heiles potential (see Fig. 4 in 
[21] for the arguments of the complex amplitude). 

Divergences of |A, C; ( lj:L )| and discontinuities of |A, c ( 3il )| of the SSPM amplitudes are seen in 
Fig. 27, too. The continuous results of the improved SPM through the bifurcations (or symmetry- 
breaking points) are also seen with the characteristic enhancement. 

Fig. 28 shows the Fourier transform of the energy-scaled quantum-mechanical level density (30). 
For a smaller a = 2.1, the diameter (2,1) orbit makes the dominant contribution to the gross-shell 
structure as the shortest POs (see the peak at r ~ 5.0). With increasing a, the amplitude of 
the circular orbit becomes again larger and shows a prominent enhancement around the bifurcation 
point (r ~ 6.2 at a^f = 7.0). Near the bifurcation, the contribution of the newborn triangular-orbit 
family (3,1), having a higher degeneracy 1C — 1, becomes also important and dominating for larger 
a > (Xuf- [The newborn (3,1) peak again cannot be distinguished from the parent circular (1,1) 
orbit near the bifurcation point a^f, like the diameter and circular orbits at a close to the HO 
limit.] 

Fig. 29 shows a nice agreement of the fine-resolved semiclassical and the quantum level densities 
Sg^e) as functions of the scaled energy e at the critical bifurcation points a = 7 and 4.25 for the 
births of the triangle-like (3, 1) and the star-like (5, 2) orbits. For the semiclassical trace formula, 
see (16), we used 

1 1 

%y, sd(e) = E = E E Ho\c) exp (-r p 2 o7 2 /4) , (32) 

K=0 K=0 po 

where 5gp^\e) is given by (29). For the quantum calculations we used the standard Strutinsky 
averaging (over the scaled energy e), finding a good plateau around the Gaussian averaging width 
7 = 2 — 3 with curvature-correction degree Ai = 6. According to Eq. (12), for the corresponding 
shell correction energy SE sc i one obtains 

6E sd /U = E ^ £ ^~m^ + n £ £ _L Sg £\e F ) . (33) 



Here, 5rii are the variations of the occupation numbers defined by the Strutinsky smoothing proce- 
dure (see, e.g., [12]). Fig. 30 shows the scaled energy shell corrections 5E sd (33) as functions of the 
particle number N 1 ^ 2 by using the standard relationship N = N(e F ), see (13), 

N = 2 g scl (e)de . (34) 
Jo 

A good plateau for the quantum calculations of the scaled energy shell corrections, see the first 
equation in (33), is realized near the same averaging parameters 7 and M. mentioned above. In 
the semiclassical calculations for the relation N(e F ), we integrated in (34) the semiclassical level 
density 

1 

07, sd (e) = 9 TF ^) + Y1 S 9\% ( £ ) > ( 35 ) 

K=0 

where g TF (s) is the scaled average part obtained within the TF approximation [4], 

OL 

9tf{£) = —7; £ , (36) 
a + 2 

and 5g\%(e) of the ISPM was determined as discussed above. For the relation N = N(e F ) we 
used specifically an averaging in small 7 = 0.1 because there is almost no sensitivity of this integral 
characteristics within the interval of smaller 7 (7 = 0.02 — 0.1). The PO sums converge for the 
averaging width 7 = 0.2 of a fine resolution of the shell structure and energy shell corrections 
with taking into account the major simplest POs [about 4 repetition numbers for the circular and 
diameter (M = 4) and a few first simplest other /C = 1 POs like (3,1), (5,2), (7,3) and (8,3)] 
which are important at a — 7. As seen from Figs. 29 and 30, we obtain a nice agreement between 
the semiclassical (ISPM, solid) and quantum (QM, dotted) results. Notice that the dominating 
contributions in these semiclassical results at the bifurcation point a = 7 are coming from the 
bifurcating circular (1, 1) and newborn (3,1) orbits, while the bifurcating circular 2(1, 1) and star- 
like (5, 2) orbits are important for a = 4.25, as expected from enhancement of their amplitudes 
seen in Fig. 27. The POs (3,1) and (5,2) yield more contributions near the bifurcation values of a 
and even more on their right in a wide region of a as mentioned above. The orbits (1,1), (3,1) [or 
(2(1,1), (2, 5)] give essential ISPM contributions of about the same order of magnitude and phase in 
the upper (or lower) panels of Figs. 29 and 30, respectively. Notice that the newborn POs (5,2), 
(7,3) and (8,3) give comparable contributions in this Figure, in nice agreement with the quantum 
Fourier spectra Fig. 28. 

The nice beating seen in these figures is explained by the interference of these orbits with the 
simple diameters of the same order in magnitude but with different phases. The latter ISPM contri- 
butions are close to the SSPM asymptotic ones near the bifurcation points a = 7 and 4.25, because 



they are far enough away from their only symmetry-breaking point at the harmonic oscillator value 
a = 2. 

4 SUMMARY AND CONCLUSIONS 

In the first part (Sect. 2) of this article we have given a short review of Strutinsky's shell-correction 
method (SCM). It was conceived as a practicable tool for numerical realizations of the theory of finite 
interacting fermion systems, developed by Migdal and collaborators, in particular at times when 
computers were not sufficiently fast to perform fully selfconsistent microscopic calculations. The 
SCM, primarily applied to atomic nuclei, combines the average binding (or deformation) energies 
from the classical liquid drop model with a microscopic shell-correction energy SE that is extracted 
from the quantum-mechanical energy spectra of realistic (deformed) nuclear shell models. It was 
and still is used world wide for very successful computations, in close agreement with experimental 
data, of nuclear masses and deformation energies, in particular fission barriers. 

In a second and larger part (Sect. 3), we have presented a semiclassical theory of quantum 
oscillations, based upon Gutzwiller's trace formula which connects the level density of a quantum 
system to a sum over periodic orbits (POs) of the corresponding classical system. This periodic 
orbit theory (POT) was extended by Strutinsky and collaborators and applied to express the shell- 
correction energy SE of a finite fermion system in terms of POs. The semiclassical trace formulae 
for SE exhibit a rapid convergence of the PO sum, due to an inverse dependence of the individual 
orbit contributions on the squares of their periods (lengths) (see Sect. 3.1). This allows one often 
to express significant features of the shell structure in terms of a few short periodic orbits. In many 
cases, the shortest POs are sufficient to describe the gross-shell features in SE. 

Bifurcations of periodic orbits can have a significant influence on a fermionic quantum system 
and leave signatures in its energy spectrum (visualized, e.g., by its Fourier transform) and hence its 
shell structure. We have (in Sect. 3.1) presented a general method to incorporate bifurcations in the 
POT, employing an improved stationary-phase method (ISPM) based on the theory of Fedoryuk 
and Maslov, and hereby overcoming the divergence of the semiclassical amplitudes of the Gutzwiller 
theory at bifurcations. The improved semiclassical amplitudes typically exhibit a clear enhancement 
near a bifurcation and on that side of it where new orbits emerge (see, e.g., Figs. 14, 15 and 27), 
which is of the order h~ l l 2 in the semiclassical parameter h. This, in turn, leads to enhanced shell 
structure effects. 

In Sect. 3.2, we have presented a semiclassical calculation of a typical actinide fission barrier us- 



ing the POT, employing a fission cavity model that uses a realistic description of the three principal 
axially-symmetric deformations (elongation, neck formation and left/right asymmetry) occurring in 
the (adiabatic) fission process (see Fig. 4). The characteristic (static) double- humped barrier (Fig. 
3) and, in particular, the sensitivity of the outer barrier to left / right-asymmetric deformations can 
be qualitatively well described by the POT (Fig. 8). The loci of minimal quantum shell-correction 
energies SE, both in particle number vs. deformation space (Fig. 5) and in two-dimensional deforma- 
tion space (Fig. 10), are correctly followed by the constant-action loci of the shortest POs. Hereby 
we observe a clear signature of period-one bifurcations of the shortest equatorial orbits which here 
were treated semiclassically using normal forms and uniform approximations. 

In Sect. 3.3, we have summarized the results of the spheroidal cavity model which had early 
been used in POT calculations for semiclassical investigations of nuclear deformations (see Fig. 2). 
This integrable system allows for analytical calculation of the PO properties. We show that the 
bifurcations of longer orbits, semiclassically treated with the ISPM, have an important influence 
on the shell structure of this system. More particularly, the bifurcated period-two and period-three 
orbits start to dominate the shell structure at medium (axis ratio rj ~ 1.5) and larger deformations. 
In the superdeformed region r\ ~ 1.6 — 2, they dominate the minima of the shell-correction energy 
5E, the more the larger the particle number (especially at r\ ~ 2), due to their enhanced semiclassical 
amplitudes, and the shortest (period-one) orbits have only a little influence (Figs. 20 and 21). The 
quantum-mechanical Fourier spectra of the spheroidal cavity exhibit a nice quantum-to-classical 
correspondence, in that enhanced Fourier signals follow exactly the lengths of the semiclassically 
enhanced classical POs (Fig. 25). This correspondence appears also in the correct description of 
the loci of large deformations in particle number vs. deformation space by the constant-action lines 
of the bifurcated period-two and -three orbits (see Fig. 24). An important reason for their strong 
enhancement at large deformations, besides the general argument given above (and explained in 
Sect. 3.1), is also the fact that the new bifurcated orbits have a larger classical degeneracy (/C = 2) 
than their parent orbits and the period-one orbits (/C = 1) (except near the bifurcations). This is 
not so in the fission cavity model (Sect. 3.2), where all POs (and all POs bifurcating from them) 
have the degeneracy K, = 1. (Another reason for differences between the two cavity models is, of 
course, the different deformation space used by them, which leads, amongst others, to a shift of 
critical deformations to larger values in the spheroidal cavity.) 

In Sect. 3.4 we have, finally, presented a class of radial power-law potentials V(r) oc r a , which 
turn out as good approximations to the popular Woods-Saxon potential in the spatial region where 
the particles are bound. The advantage of the power-law potential is that, in spite of its diffuse 



surface, the classical dynamics scale with simple powers of the energy, which makes it particularly 
easy for POT calculations. The quantum Fourier spectra yield directly the lengths of the leading 
classical POs (Fig. 28). We have developed semiclassical trace formulae for this potential and 
study various limits of the power a (harmonic oscillator potential for a — 2, cavity potential for 
a — > oo). Bifurcations are, again, treated in the ISPM leading to semiclassical enhancement of 
the orbit amplitudes. The trace formulae are shown numerically to give good agreement with the 
quantum-mechanical level density oscillations (Fig. 29) and shell-correction energies (Fig. 30). 

In conclusion, we may state that the semiclassical theory (POT) is well capable of explaining the 
main features of quantum shell structure in terms of a few classical periodic orbits. Bifurcations of 
POs are not simply an obstacle of the semiclassical theory, but they leave clear signatures both in 
the quantum Fourier spectra and the locations of minima of the shell-correction energy SE plotted 
versus particle number and deformations. 

In our future work we intend to further study finer shell structures due to longer orbits, applying 
both the ISPM and uniform approximations to treat their bifurcations. In particular, we want 
to understand the "quantization" into isolated minima that occurs along the valleys of both the 
ground-state and the secondary (isomer) minima and, in particular, we want to better understand 
the superdeformed shell structure. In the region of super deformation in the spheroidal cavity, we 
see clearly from Fig. 24 that this is due to the bifurcated period-two and -three orbits. It will be 
interesting to study the analogue situation in the fission cavity model by investigating the higher- 
period bifurcations occurring there. One complication in that model is that we have to vary three 
deformation parameters (c, h, a), instead of only one (rj) for the spheroid, in order to fully cover the 
fission barrier landscape. This will naturally lead to a much larger variety of possible bifurcations. 
Another object of our studies will be to understand semiclassically the transition regions between 
the downgoing ground-state valleys and the initially upgoing secondary-minimum valleys (see Figs. 
5 and 23). They are well described separately by the shortest meridional and equatorial orbits, 
respectively, but in the transition regions (1.2 < c < 1.35 and 1.2 < r\ < 1.5, respectively) we 
expect further bifurcations and corresponding new orbits to play a role. As a more technical 
aspect, we also want to perform detailed comparisons between the results obtained with the ISPM 
and the uniform approximations for specific bifurcations. 

On a longer time scale, we intend also to include pairing and spin-orbit interactions into the 
POT calculations, in order to come closer to a realistic description of nuclei, and to extend the POT 
towards dynamics (e.g., studying inertial moments, friction coefficients, etc.). We keep in mind and 
emphasize, however, that the POT is not intended as a substitute for a fully quantum-mechanical 



theory, but rather a tool for a better understanding of its results in terms of classical pictures. 



Acknowledgement 

We are very grateful to our late teachers, Profs. A. B. Migdal and V. M. Strutinsky, for so many 
creative ideas and collaborations. We also thank Profs. K. Matsuyanagi and E. E. Saperstein for 
fruitful collaborations and many useful discussions. 



Appendix A: 
2D RADIAL POWER-LAW POTENTIAL 



Al. Classical Dynamics in Generic 2D Circular Hamiltonians 

For an isotropic 2D potential V(r), the Hamiltonian H in the canonical phase-space variables 
{r, <p,p r ,P v } writes 

ff = i(^ + S) + v >» = £ - < A1 > 

Here, p v = £ is the angular momentum, p> is the cyclic variable, 

ip = u v t + p>' , (A2) 

tp' — ip(t — t' — 0) is the initial condition, 



£2 £2 

p r (r) = Y p 2 (r) = y 2m [E — V(r)\ = mr, p(r) = y 2m [E — V(r)\ . (A3) 

In particular, we shall consider the power-law potential (23) from [28]. 

Integrating the first equation in (A3) one has also the radial trajectory r(t), 

m ^ dp = t , (A4) 

Jr> ^2m [E - V(p)] - £ 2 /p 2 

where r' — r(t — t' — 0) as the initial condition for the radial coordinate r. For the radial action I r 

in polar action-angle variables over one period (t = T) with p' = <p" (r ; = r") one has 

If 1 r Tmax I f2 

I r = —j )Pr dr = - drJ2m[E-V(r)]-- = I r (EJ), 

I<p = j P v dp = £ , (A5) 

where r min and r max are the turning points for a given energy E, which are solutions of the equation: 

I 2 

T(r, £) = 2m [E — V{r)\ - — = . (A6) 

For frequencies one has 

_ dH 1 (m f rmax dr \ 

Ur= dT r = dI r /dE = J rmin ^2m [E - V{r)]^Wp ) ' 

where I r = I r (E,£) is the surface energy (A5). Thus, the periodic-orbit equation writes 

m -^- dJ r(E,e) _e r— dr = n^ 

u r ~ dt ~ 7T J rmm r y2m [E — V(r)] — £ 2 /r 2 n r ' 

n v and n r are co-primitive integers. The energy surface I r = I r (E,£) is simplified to function of 
one variable I. The solutions to the PO equation (A8) for i = t{n r ,n v ) define the JC = 1 families 
of orbits M(n r ,n v ), see [15]. 



A2. Circular Orbit in Power-Law Potential 



In the radial power-law potential (23), the isolated circular orbits cause successive bifurcations with 
increasing a. At the energy E, the radius of a circular orbit is given by 



P ■ 2E 

r °~ \(2 + a)U 



l/a 



4 = P{rc)r c , T c = y~ < 0, 



T c is defined by (A6). The radial frequency Q c = ou r (r c ) is given by 

2E 



a 



a(2 + a)U 
mR 2 



(2 + a)U 



> 0. 



From (A7) at £ = 4 by using (A9) for r c and 4 one finds [28] 

2E ^ l/2-l/a 



u c = u (p (£ = 4) 



From (A10) and (All), one has the equation, 



all 



mR 2 \(2 + a)U 



given explicitly by [28] 



V2 + a 



n r 



(A9) 



(A10) 



(All) 



(A12) 



(A13) 



Thus, for the bifurcation diffuseness parameter a bi j- where the one-parametric orbits (n r ,n v ) are 
born one arrives at the explicit expression: 



«w/ = n ll n l ~ 2 



(A14) 



A3. Trace formula for 2D Circular Power-Law Potential Model 

From the ISPM procedure presented in sect. 3.1, we obtain contributions of the /C = 1 families and 
isolated circular orbits into the semiclassical level density. For a K, — 1 family, we choose the angular 
momentum t and angle <p conjugate to £ for integration variables in the phase-space trace formula 
(3). The integration over a cyclic angle (p generates only the factor 2ir. The remaining integration 
over £ is carried out within the simplest ISPM by expanding a phase function $ in powers of £ — £* 
up to quadratic terms. [The expansion of pre-exponent factor of integrand in powers of £ — £* is cut 
at zero order in this approximation.] 

Thus, we obtain the amplitude Ap~J of Eq. (26) as 



with the period 



T po = ^ = ^ , (A16) 

0J r Up 



K po is the curvature of the energy surface, 

K _dH r _ dHAEj) df(£) 

d/J ~ dP ~ ~dT' [An) 

f(£) is the ratio of frequencies defined in the first three equations of (A8). The arguments of the 
generalized complex error function in (A15), see also (10), is given by 



Z, 



= ^-inMn r \K po \/hb po {£± - i po ) , (A18) 



where b po = 1 for all /C = 1 polygon-like PO families, except for the diameters for which one has 
(see [21]) 



b M (2,i) = 1 - 2 ex P 



a — 2 



2A 



D 



^ixMn r K D jti 



(A19) 



Kb is the diameter curvature. For simplicity, the finite integration interval of the angular momenta 
was split into two parts, — £ c < £ < and < £ < £ c , where £ c is the angular momentum of 
a circular orbit as mentioned above. There are the symmetric stationary points ±\£*\ related to 
the clockwise and anticlockwise motion of the particle along the PO in these two parts of the 
phase space. They give equivalent contributions to the amplitude owing to the independence of the 
Hamiltonian of time. Thus, we have reduced the integration region to < £ < £ c accounting for this 
time- reversibility symmetry simply by a factor 2 in (A15). Therefore, for the polygon-like (n r > 2) 
POs, one has £=0 and £ + = £ c . In the HO limit a — > 2, only the diameter (2,1) and the circle (1,1) 
(with repetitions) survive and form /C = 2 families in the HO potential. It is natural to assume that 
the whole phase-space in the HO limit can be split into the two equivalent parts, < £ < £ c /2 for 
the diameter and £ c /2 < £ < £ c for the circular orbits, for the angular momenta and correspondingly 
for the conjugate angles. Thus, we arrive at the integration limits (A18) with (A19). The same 
prescription has been successfully applied to the integrable version of the Henon-Heiles potential for 
contributions of orbits B and A in the same HO limit [21]. For the Maslov index of the considered 
/C = 1 family of POs, one obtains 

<7« = 2Mn r , (A20) 

and for the constant phase independent of the orbits one has (fid = — 7r/4. 

Extension of the limit £± — > ±oo exactly reproduces the asymptotic Berry- Tabor formula [16, 4] 



4(1) - d po T po / A91 x 

p ° nhVyMnr \K P ' ■ 



po\ 



where d po takes into account the discrete degeneracy, d po = 1 for diameters M(2, 1) (n r = 2n v ) and 
2 for all other POs (n r > 1n^) [4]. In the circular billiard limit (a — > oo), the action S po (E) is given 
by 

S po {E) = pL po , p = V2^E , (A22) 
and the curvature K po from (A17) is obtained explicitly by 



1 



K = 
po ixpR sin ■ 



(A23) 



Substituting these quantities, Eqs. (A22) and (A23), into (A21), from (26) and (A20) one obtains 
the Balian-Bloch trace formula for the circular billiard [4]. 

For the contribution of circular orbits, one should perform the two exact phase-space integrations 
over radial variables within ISPM. Finally, we arrive at the result: 



A {0) 



eri(Z^,Z^)eri(Z^,Z 



rc ' rc / ' 



2nhy/F CtM 

where T c is the period of a primitive circular orbit. F C>M is its stability factor [13] given by 



(A24) 



c,M 



4 sin 



irMy/a + 2 



(A25) 



The arguments of the complex error functions in (A24) are given by 



= 

p c 



Z (±) 



--7i M Va + 2 
n 



i F, 



c,M 



where K c is the circular-orbit curvature, 



K, 



(£ ± - 4) , h = 4, l- = 4/2 , 

(±) , 0(+) = vr, 0( r ) = o , 



(A26) 



(a + l)(a-2) 



K c = _ _ 

12(V^T2) 3 4 



(A27) 



The time-reversibility of the Hamiltonian was similarly taken into account as explained above. For 
the Maslov index one has a p °J = AM and 4> d = 0. (The detailed derivation will appear elsewhere.) 
In the asymptotic limit of the finite non-zero integration boundaries — > —oo and 0+ — > oo, i.e. 
far from any bifurcations a bi j (also from the harmonic oscillator symmetry breaking at a = 2), the 
expression (A24) tends to the amplitude of the Gutzwiller trace formula for isolated orbits [13, 4]. 
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Figure 1: Schematic picture showing the relation of local level density at the Fermi energy to the 
binding of the system. Left side: Low level density - more bound; right side: high level density - 
less bound. 




Figure 2: Neutron energy shell-correction SE n (N,rj) for a realistic nuclear Woods-Saxon potential 
(with spin-orbit term), versus neutron number N and axis ratio r\ of the spheroidally deformed 
potential (equidistance: 2.5 MeV, areas with negative values are shaded). Dots indicate experimen- 
tal deformations. The heavy bars are the semiclassical predictions for the loci of the ground-state 
minima using the leading classical periodic orbit families in a spheroidal cavity. 




Figure 3: Left: Schematic double- humped fission barrier of a typical actinide nucleus. Note the 
lowering of the outer barrier due to left-right asymmetric shapes. Right: Maximum probability 
amplitudes (schematic) of the two leading s.p. states responsible for the asymmetry effect. 
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Figure 4: Axially symmetric nuclear shapes in the (c,h,a) parameterization of [12]. Dashed lines 
for a^O; the sequence with h = a = corresponds to the shapes obtained in the LDM [33]. 
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Figure 5: Plot of quantum-mechanical shell-correction energy SE versus cube-root of particle num- 
ber, iV 1//3 , and elongation c (along h = a = 0) in the cavity model. The contours lines are for 
constant values of SE (white: positive values, gray to black: negative values. The heavy lines in- 
dicate the loci of constant actions of the leading POs. Dashed- dotted lines: meridional triangular 
orbits (3,l,l)s; narrow lines: diameter orbits in equatorial planes (2,1)EQ (solid) and in parallel 
perpendicular planes (2,1)AQ (dashed); broad lines: triangular orbits in equatorial planes (3,1)EQ 
(solid) and in parallel perpendicular planes (3,1)AQ (dashed). The horizontal dotted line at iV ~ 180 
corresponds to the situation with the isomer minimum at the correct deformation [12] c ~ 1.42 of 
the real nucleus 240 Pu. 




Figure 6: Fourier spectra of the fission cavity model with h = a = for five values of c: amplitude 
of Fourier transform of the quantum spectrum versus length L (in units of Rq) of the classical 
POs. Short arrows: POs lying in planes orthogonal to symmetry axis; long arrows: POs lying in 
meridional planes containing the symmetry axis (labels as in Fig. 5; see text for more details). 
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Figure 7: Quantum-mechanical (solid line) and semiclassical (dashed line) level density Sg versus 
elongation c (along h = a = 0) in the fission cavity model taken at the Fermi wave number 
kp = 12.1/ Rq. (Gaussian averaging over wave number k with width 7 = 0.6/i? -) 




Figure 8: Perspective view of the semiclassical outer fission barrier versus elongation c and left- 
right asymmetry a for h — 0. To the left, the shapes corresponding to the points A, B and C in 
the deformation energy surface are displayed; the vertical solid (dashed) lines indicate the planes 
containing the stable (unstable) POs. 




Figure 9: Squared amplitudes of the Fourier transforms of Sg(k) obtained (as in Fig. 7) quantum 
mechanically (solid line) and semiclassically (dashed line) for c=1.5, h—0, a— 0.12. The insert 
shows the cavity shape and the location (vertical line) of the plane containing the shortest POs 
with reflection number p. The abscissa shows the length L of the POs in units of the spherical- 
cavity radius R ; see the text for further details. 




Figure 10: Contour plots of the shell-correction energy 5E versus c and a. Upper panels: for h — 0, 
lower panels: for h = —0.075. Left panels: results of quantum-mechanical SCM calculations with 
realistic nuclear shell model potentials [12] (shown is the shell correction of the neutrons); right 
panels: semiclassical POT results with the fission cavity model described above. 



Figure 11: POs M(n u ,n v ,n v ) in the spheroidal cavity model: 2D orbits. 




Figure 12: POs M(n u ,n„,n ¥ ,) in the spheroidal cavity model: EQ orbits. 




Figure 13: POs M(n u ,n v ,n v ) in the spheroidal cavity model: 3D orbits. 




Figure 14: Amplitudes ^.5,2,1 ■ Solid and dotted lines show the quantum and semiclassical ISPM 
results, respectively. 




Figure 15: Amplitudes Aq^i- Solid and dotted lines show the quantum and semiclassical ISPM 
results, respectively. 
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Figure 16: Oscillating part of level density 5g(k) in the spheroidal cavity (in units of 2mR 2 /h 2 , 
where R is the radius of the cavity at rj — 1), plotted versus kR, for three specific deformations r\. 
The Gaussian averaging parameter is 7 = 0.3/ R. 
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Figure 17: The same as in Fig. 16 but for three deformations in the super deformed region. 
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Figure 18: Shell-correction energy SE (in units of h 2 /2mR 2 ) vs. cube root of particle number iV 1//3 
for the same three deformations as in Fig. 16. 
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Figure 19: The same as in Fig. 18 but for the three deformations of Fig. 17. 
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Figure 20: The same as in Fig. 19 in the super deformation region, but the semiclassical results 
(ISPM) here include only the bifurcated period-two and -three orbits specified near the dotted lines 
"ISPM". 




Figure 21: Shell-correction energy 8E (in units of H 2 /2mR 2 ) vs. cube root of particle number iV 1//3 
for the characteristic deformations rj = 1.5, 1.732... and 2.0. Solid lines are the quantum results, 
dots represent the ISPM results including the bifurcating period-two orbits. The thick solid lines 
show the ISPM results obtained using only the shortest EQ and meridional orbits (see text for 
discussion). 




Figure 22: Fourier transform \F(L)\ (14) of the quantum level density g(k) of the spheroidal cavity, 
plotted vs. length variable L at the critical deformations rj = 1.5, 1.62, 1.73 and 2.0. It shows a 
correspondence to the lengths of classical periodic orbits and clear dominance of longer bifurcated 
period-two POs above the shortest period-one POs: arrows near the peaks mark the lengths of the 
most important POs; some peaks at the bifurcation deformations show the sum of both newborn 
and parent bifurcating orbits like in Fig. 6. 
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Figure 23: Contour map of shell correction energy SE of the spheroidal cavity versus deformation rj 
and cube- root of particle number N 1 ^ 3 . Thick and thin solid lines represent constant-action curves 
of the classical shortest meridian POs; thick and thin dashed lines show the shortest EQ POs; all 
are specified on right and are dominating the periodic-orbit sum for the gross-shell structure. 
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Figure 24: The same as in Fig. 23 but with constant-action lines of the bifurcating POs. Three kinds 
of curves (dotted, dashed and solid) are shown: Normal and heavy dotted lines present the constant- 
action curves of the bifurcating 2(2,1)EQ and newborn (4,2,1)2DH orbits; similar dashed lines the 
(5,2)EQ and (5,2,1)3D orbits, and solid lines show the 2(3,1) and (6,2,1)3D orbits, respectively. All 
newborn 2DH and 3D orbits dominate the fine shell structure at large enough deformations. 




Figure 25: The Fourier transform |-F(L)| (14) of the quantum level density g{k) is plotted in 
grayscale. It shows a clear correspondence to the lengths of classical periodic orbits: dotted curves 
represent EQ orbits, and dashed curves stand for meridional orbits with elliptic caustics 2DE. Solid 
curves are either meridional orbits with hyperbolic caustics 2DH or 3D orbits, which are generated 
via period-multiplying bifurcations of EQ orbits. The black dots show the bifurcation points. 
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Figure 26: Some short periodic orbits in the power-law potential (23) for several values of a. 




Figure 27: The scaled amplitudes |^4 SCiP o| as functions of a for the once repeated (M = 1) circular 
[(1, 1), see (30)], one-repeated (M = 1) diameter (2, 1) and triangle PO (3, 1), in units of e 1 / 2 at the 
scaled energy e = 10. Dashed thin (circular (1,1)SSP, diameter (2,1)SSP) and thick (3,1)SSP lines 
are SSPM results, and solid (thick circular (1,1)ISP and thin diameter (2,1)ISP) as well as dotted 
(3,1)ISP lines are ISPM results. 
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Figure 28: Moduli of the Fourier transform |F(r)| of the energy-scaled quantum level density (30), 
plotted for several values of a. 




Figure 29: Oscillating part of scaled level density versus scaled energy e at the bifurcation points 
a = 7 and 4.25 in units of e 1 ^ 2 . QM (solid line) is the quantum- mechanical result (using Strutinsky 
averaging with parameters 7 = 3, M = 6); ISPM (dotted line) is the semiclassical result. In both 
(dimensionless) width 7 = 0.2 was used in the Gaussian coarse-graining over the scaled 

energies. 




Figure 30: Shell-correction energy 5E vs. particle number parameter iV 1//2 at a = 7 and 4.25 in 
units e 3 ( a ~ 2 )/[ 2 ( a + 2 )l . QM (solid line) is the quantum-mechanical result using the shell-correction 
method (Strutinsky averaging as in Fig. 29 above); ISPM (dotted line) is the semiclassical result 
using a small coarse-graining width 7 = 0.1 in the calculation of N(ef) by (34) and (25). 
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